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1.  INTRODUCTION 


The  main  objective  of  this  project  is  to  develop  a  computer-aided  multivariable 
control  system  design  technique  by  matching  the  frequency  responses  of  the  compensated 
closed  loop  system  and  a  given  ‘desired’  transfer  function  matrix  over  a  frequency  interval 
of  interest.  The  motivation  for  this  approach  is  to  come  up  ‘quickly’  with  a  candidate  for 
the  compensation  which  could  be  ‘fine  tuned’  for  practical  applications. 

The  proposed  MIMO  (multi— input  multi-output)  control  design  technique  is  a 
generalization  of  the  SISO  (single-input  single-output)  technique  of  the  Principal 
Investigator  [13].  An  aircraft  application  of  the  SISO  technique  was  presented  by  Gentry 
and  Pujara  in  [4], 

Under  the  MIMO  design  technique  presented  in  this  report,  one  starts  with  a  trans¬ 
fer  function  matrix  of  a  stable  plant  and  the  objective  is  to  design  a  MIMO  controller  so 
that  the  overall  closed  loop  transfer  function  matrix  is  ‘close’  to  a  given  desired  transfer 
function  matrix.  The  desired  transfer  function  matrix  is  synthesized  to  incorporate  a  given 
set  of  specifications.  As  was  mentioned  earlier,  the  plant  has  to  be  stable  and  if  it  is 
unstable,  then  it  has  to  be  stabilized  before  applying  the  proposed  technique. 

Under  the  proposed  technique,  the  controller  parameters  are  obtained  by  minimiz¬ 
ing  a  weighted  mean  square  error  between  the  frequency  responses  of  the  compensated 
closed  loop  system  and  a  desired  closed  loop  system.  The  weighting  function  in  the  error 
mentioned  above  is  chosen  in  such  a  way  that  the  accuracy  of  the  error  is  maintained  as  far 
as  possible  and  in  addition  to  this,  the  parameters  of  the  controller  turn  out  to  be  solutions 
of  linear  algebraic  equations.  This  last  aspect  of  the  algorithm  is  extremely  important  in  as 
much  as  this  leads  to  a  simple  computer-aided  package  to  implement  the  algorithm. 

The  program  to  implement  the  algorithm  has  been  written  in  Fortran  77  and  is 
quite  simple  to  use.  All  that  the  user  has  to  do  is  to  input  the  coefficients  of  the  plant 
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transfer  function  matrix,  the  limits  of  frequency  interval  and  the  order  of  the  controller. 
These  frequency  intervals  are  identified  by  plotting  the  magnitude  response  of  the 
components  of  the  desired  transfer  function  matrix.  As  output,  the  user  gets  the 
parameters  of  the  controller,  the  frequency  responses  of  the  compensated  closed  loop 
system  and  the  desired  transfer  function  matrix.  If  this  is  satisfactory,  one  is  done.  If  not, 
then  the  user  increases  the  order  of  the  controller  and/or  ‘tinkers’  with  the  limits  of 
integration  for  the  error  function  until  satisfactory  results  are  obtained.  This  ‘trial  and 
error’  is  not  excessively  time  consuming  because  it  is  a  computer-aided  technique.  Four 
examples  are  solved  in  this  report  for  illustration  with  excellent  results.  Example  1  is  the 
same  example  studied  by  Chen  [1]  and  the  comparison  of  the  results  by  the  proposed 
technique  with  those  of  Chen  [1]  shows  that  the  proposed  technique  gives  letter  results. 
The  three  parts  of  the  second  example  are  related  to  YF— 16  CCV. 

In  addition  to  this,  a  modest  ‘robustness’  analysis  of  the  design  for  the  YF— 16  CCV 
control  system  design  obtained  by  the  proposed  technique  in  Example  2  has  been  carried 
out.  Only  one  parameter  in  the  plant  after  it  was  stabilized  was  varied.  By  using  the  same 
controller  as  obtained  for  the  nominal  plant,  it  was  found  that  for  a  larger  range  of  values 
for  this  parameter,  the  inner  loop  maintains  its  stability  and  the  overall  closed  loop  system 
also  maintains  an  excellent  match  with  the  desired  closed  loop  transfer  function  matrix. 
This  has  been  quantified  by  drawing  a  graph  of  the  error  function  (between  the  frequency 
responses  of  the  compensated  system  and  the  desired  system  over  the  frequency  interval 
[0,10]  rad.)  against  the  changes  in  the  parameter.  The  frequency  interval  [0,10]  rad/sec 
was  picked  as  this  interval  is  of  great  significance  from  the  point  of  aircraft  flying  qualities. 
In  addition  to  this,  some  sample  graphs  (over  the  range  of  the  parameter)  for  magnitude, 
phase  and  time  response  have  been  given  which  portray  good  match  and  reinforce  the 
results  predicted  by  the  error  function.  Evidently,  a  lot  more  time  and  effort  is  required  to 
do  a  complete  'robustness  analysis  for  the  example  under  consideration  here.  But  this  is  a 
reasonably  good  start. 
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It  should  be  pointed  out  that  the  design  technique  proposed  in  this  report  is  only  for 
nominal  plant  and  no  disturbance-rejection  or  noise— suppression  considerations  were 
included  in  the  error  criterion.  It  is  a  complicated  matter  and  the  inclusion  of  these 
concerns  in  the  design  methodology  is  left  for  future  research  work. 

A  program  written  in  Fortran  77  in  the  form  of  a  floppy  diskette  has  already  been 
submitted  to  Mr.  Tom  Gentry,  the  Project  Engineer. 
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2.  A  BRIEF  REVIEW  OF  THE  FREQUENCY  MATCHING  TECHNIQUE  FOR 
SINGLE-INPUT  SINGLE-OUTPUT  CONTROL  SYSTEMS 


2.1  Introduction 

The  objective  of  this  chapter  is  to  review  the  frequency  matching  technique  for  the 
SISO  (single— input  single-output)  control  system.  This  technique  uses  the  idea  of  a 
weighting  function  to  ‘linearize’  the  problem,  which  has  been  in  the  literature  for  quite 
some  time.  This  idea  wras  originally  used  by  Levy  [10]  for  complex  curve  fitting,  and  then 
used  by  several  researchers  [14,15,16,17,18,19]  for  model  reduction  problems  and  digitizing 
continuous  systems.  In  addition  to  this  linearizing  idea,  we  have  added  several  other 
features  in  the  error  criterion  to  make  it  a  viable  and  practical  design  tool.  Under  this 
technique,  the  controller  parameters  are  obtained  by  matching  the  frequency  responses  of 
the  compensated  closed  loop  system  and  a  ‘desired’  closed  loop  system.  The  ‘desired’ 
transfer  function  is  synthesized  to  meet  a  given  set  of  specifications,  which  for  example 
could  include  flying  qualities  criteria  for  an  aircraft  control  system  design. 

The  error  is  defined  in  such  a  way  that  the  controller  parameters  turn  out  to  be 
solutions  of  linear  algebraic  equations.  This  is  quite  significant  [2,6]  when  contrasted  with 
the  results  where  the  controller  parameters  turn  out  to  be  solutions  of  nonlinear  algebraic 
equations  which,  in  general,  are  harder  to  solve. 
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2.2  The  SISO  Frequency  Matching  Technique 


KO 


Figure  2.1:  Block  diagram  for  SISO  control  systems 


Consider  a  typical  unit  feedback  control  system  as  represented  by  its  block  diagram 
as  shown  in  Figure  2.1  where  G(s)  is  the  plant,  assumed  stable,  and  C(s)  is  a  controller.  It 
is  desired  to  synthesize  a  controller  C(s)  so  that  the  compensated  closed  loop  transfer 
function  is  ‘close’  to  a  certain  ‘desired’  closed  loop  transfer  function  over  a  frequency 
interval  of  interest.  Suppose  D(s)  is  a  ‘desired’  closed  loop  transfer  function  which  has 
been  synthesized  to  satisfy  a  given  set  of  specifications.  Suppose 


c(i)  =  - ■  m£n 

i+5>- 


where  a’s  and  b’s  are  unknown  real  parameters.  Define 


G(;'u?) 


L\(y>)  +  jLj(u>) 
Mi  (<*0  +  jMa{v) 


(2.2) 


D(JU>'  Pi(uf)  +  jP,W 


(2.3) 


where  L^w),  L2(w),  M^w),  M2(u/),  Ni(cj),  N2(u/),  P](u>)  and  P2(w)  are  polynomials  in  u. 
Then  the  frequency  response  of  the  compensated  closed  loop  transfer  function  can  be 
written  as 
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»=0 _ 

>+ 


L\(ui\  4  jL2(u>) 

,  4  jMi(u>) 


F{ju>)  = 


£  a,(ju>y 

1  ,  I  .=0 _  Fj(u>)  -f  jL2( u>) 

1+ 


As  in  [13],  the  weighted  mean  square  error  between  the  frequency  responses  of  F(jw)  and 
D(ju>)  can  be  written  as 

E  =  I  J[A(**>)  +  jP*{u)]\Li{w)  +  jlj(a>)]  ^a.O'u;)' 

J  Ui\  .  _ 

1  l»=0 

-[7V,(u;)  +  jA^2(u))3 /[A/^w)  +  jA/2(u))]  1  +  ^2  htiu,y  (2.5) 
t  1=0 

+  [Z,i(w)  +  <7'Ia(ta>)]  f^a,(iw),j  ||  du> 


where  [uq,^]  is  the  frequency  interval  of  interest. 

As  in  [13],  the  minimization  of  E  with  respect  to  the  controller  parameters,  results 
in  linear  algebraic  equations  involving  a’s  and  b’S.  These  equations  are  written  as  a  system 
of  equations  in  matrix  form  as  was  described  in  [13].  The  complete  derivation  and 
application  of  this  frequency  matching  technique  for  the  SISO  control  systems  algorithm 
can  be  found  in  [13]  and  [4],  respectively. 


2.3  Summary 

In  this  chapter,  a  computer-aided  technique  of  designing  a  SISO  control  system  is 
reviewed.  The  controller— parameters  are  obtained  by  matching  the  frequency  responses  of 
the  compensated  closed  loop  system  and  a  ‘desired’  closed  loop  system  over  a  certain 
frequency  interval  of  interest.  The  technique  is  easy  to  simulate  as  the  controller  para¬ 
meters  turn  out  to  be  solutions  of  linear  algebraic  equations.  A  much  more  detailed 
explanation  of  the  SISO  frequency  matching  technique  can  be  found  in  [13].  In  addition  to 
this,  [4,13]  is  another  good  reference  as  an  illustration  of  potential  applications  of  the 
technique  for  aircraft  control  systems  design. 
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3.  MULTIVARIABLE  FREQUENCY  MATCHING  TECHNIQUE 


The  objective  of  this  chapter  is  to  generalize  the  frequency  matching  technique  for 
the  SISO  control  system  design  case  to  the  MIMO  case.  We  begin  with  a  brief  intro¬ 
duction  to  the  multivariable  function  matching  technique  algorithm.  Then  we  do  the 
problem  formulation  and  then  give  the  complete  derivation  of  the  algorithm.  Finally,  we 
give  some  explanatory  remarks  regarding  the  application  of  the  technique. 

3.1  Introduction 

The  main  thrust  of  this  report  is  to  develop  a  computer-aided  control  system  design 
technique  for  linear  multivariable  continuous  control  systems  via  frequency— matching.  As 
in  the  SISO  case,  this  is  accomplished  by  matching  the  frequency  responses  between  the 
compensated  closed  loop  system  and  a  ‘desired’  closed  loop  system  by  means  of  a  weighted 
mean  square  error.  A  suitable  generalization  of  the  error  function  for  the  MIMO  case 
would  be  defined  to  obtain  the  parameters  of  the  controller. 

Wc  propose  to  obtain  the  controller  parameters  by  minimizing  a  weighted  mean 
square  error  between  the  frequency  responses  of  the  compensated  closed  loop  transfer 
function  matrix  and  the  ‘desired’  closed  loop  transfer  function  matrix  over  a  certain 
frequency  interval  of  interest.  The  ‘desired’  closed  loop  transfer  function  matrix  is 
synthesized  to  meet  a  given  set  of  specifications.  A  similar  idea  of  using  a  ‘desired’  transfer 
function  matrix  for  multivariable  control  system  design  has  been  used  by  Lehtomaki,  Stein 
and  Walls  [9],  Chen  [1],  etc.  As  in  the  SISO  case,  the  error  is  defined  in  such  a  way  that 
the  controller  parameters  turn  out  to  be  solutions  of  linear  algebraic  equations.  To  make 
the  technique  fairly  general,  intervals  of  frequency  for  different  elements  of  the  transfer 
function  matrix  could  be  used.  The  derivations  were  quite  involved  and  intricate,  but 
eventually  quite  useful. 
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3.2  Problem  Formulation 


Consider  a  typical  multivariable  feedback  control  system  of  order  n  represented  by 
the  following  block  diagram 


Figure  3.1:  Multivariable  feedback  control  system 

where 

r(s)  =  vector  of  reference  inputs,  of  order  n. 
y(s)  =  vector  of  outputs,  of  order  n. 

G(s)  =  n  x  n  transfer  function  matrix  of  the  plant. 

H(s)  =  n  x  n  transfer  function  matrix  of  the  feedback  element. 

C(s)  =  n  x  n  transfer  function  matrix  of  the  controller 

which  has  to  be  synthesized  to  meet  the  design  specifications. 
The  closed  loop  transfer  function  matrix  F(s)  for  this  system  is  given  by 


F(«)  =  (I.  +  G(.)C(.)H(*))-'  G(j)C(.) 


(3.1) 

where  In  is  the  identity  matrix  of  order  n.  Suppose  D(s)  is  a  ‘desired’  n  x  n  closed  loop 
transfer  function  matrix  which  has  been  synthesized  to  satisfy  a  given  set  of  specifications. 
These,  for  example,  could  include  flying  qualities  criteria  for  aerospace  application. 
Suppose 


G(-)  = 

C(-)  = 


[*»(«)) 
M  ’ 

M«)l 

a(s)  ’ 


"w  -  ‘t#  • 


D(a)  -  d$f  ’ 


(3.2) 

(3.3) 
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p(«,  k)  <  9  > 


and 


p(*>) 

CiJs)  -  £ 


a<l>V 


a(3)  =  1  +  5  ,  ■®F 


r  —  1 


1  <  : ,  k  <  n 


(3.4) 


where  g’s,  h’s,  c’s.  d’s  and  a(s)  are  polynomials  in  s  with  real  coefficients.  The  controller 
C(s)  has  a  common  denominator  a(s),  with  and  bz,  for  0  <  r  <  p(i,k)  and  1  <  i,  k  <  n, 
and  1  <  z  <  q,  respectively,  as  the  real  unknown  parameters.  In  this  chapter,  a  notation  of 
the  form  [gik(s)]  will  stand  for  annxn  matrix  with  gik(s)  as  its  (i,k)th  entry. 

Now  we  define  a  mean  square  error  between  the  frequency  responses  F(jw)  of  the 
compensated  closed  loop  system  and  D(jw)  of  the  ‘desired’  closed  loop  system,  to  determine 
the  unknown  parameters  in  the  controller.  This  is  defined  as  follows 


£(.,t)  =  EE  if.*Ow)  -  a.M 


du) 


(3.5) 


where  Fik(jw)  and  Dik(jw)  are  the  (i,k)th  elements  of  the  matrices  F(jw)  and  D(jw), 

/ 

respectively.  Also  [wik,wik]  is  the  frequency  interval  of  interest  for  the  (i ,k)th  element  in 
the  ‘desired’  transfer  function  matrix  D(s)  which  usually  is  the  bandwidth  of  that  partic¬ 
ular  transfer  function.  Also  a  and  b  are  the  vectors  of  the  unknown  controller  parameters, 
a!1,k)  for  0  <  r  <  p(i,k),  1  <  i,k  <  n,  and  bz  for  1  <  z  <  q,  respectively.  An  optimal  set  of 
{a[1,k',bz}  for  0  <  r  <  p(i,k),  1  <  i,  k  <  n,  and  1  <  z  <  q,  can  be  obtained  by  minimizing 
E(a,b). 

It  can  easily  be  checked  that  Fik(jw)  is  the  (i,k)th  element  of 

g{]*')h(ju>)a{ju>){6lk  jg,t(ju;):!clk(ju>)][/ta(ju>)] 

g{ju;)h{ju>)a{ju>) 

where  ^ik  is  the  Kronecker  delta.  Now  to  obtain  the  controller  parameters,  we  minimize 
the  error  function  E(a,b)  as  in  equation  (3.5)  by  equating  to  zero  the  partial  derivatives  of 
E(a,b)  with  respect  to  and  bz,  for  0  <  r  <  p(i,k),  1  <  i,  k  <  n,  and  1  <  z  <  q.  This  way 
one  will  get  a  set  of  complicated  nonlinear  algebraic  equations. 


[9*k(jv)][c,k(j<*>)} 

g(ju;)a{ju;) 


(3.6) 
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3.3  The  Algorithm 

We  will  now  show  that  if  the  error  E(a,b)  is  suitably  modified,  while  retaining  a 
reasonable  degree  of  accuracy  of  approximation,  the  multivariable  frequency— matching 
technique  will  lead  only  to  solving  of  much  simpler  linear  algebraic  equations  involving  the 
controller  parameters.  We  feel  that  this  is  a  key  point  in  the  development  that  follows. 
We  want  the  frequency  response  of  the  compensated  closed  loop  transfer  function  matrix, 
F(ju>)  to  match  the  frequency  response  of  the  ‘desired’  closed  loop  transfer  function  matrix, 
D(ju/)  over  a  frequency  interval  of  interest.  In  symbols,  we  wish  to  Achieve 


[I„  +  G(jtv)C(jw)H(ju))]-1  G(ju>)C(ju>)  ~  D(jw) 
G(ju>)C(ju>)  a  (I„  +  G(;u»)C(ju>)H(ju>)]  D(ju/) 


for  all  u)  in  the  frequency  interval  of  interest. 

Making  appropriate  substitution,  we  thus  wish  to  achieve 


|g»*(j^)][Cil.(ju>)] 

\g(jv)h(ju)a{ju>)[6,k}  +  lg,fc(j^)Mj*))lMM 
*  [  ”  g{ju>)h(j»)a{ju>) 


d(ju>) 


(3.7) 

(3.8) 


(3.9) 


Now  we  cross  multiply  in  the  above  Equation  (3.9).  Thus  we  see  that  we  wish  to  achieve 

*  \g{j*)h(j*)a{ju>)[6ik}  +  (**(>")) 

(6AW) 

for  all  ui  in  the  frequency  interval  of  interest.  This  has,  in  some  sense,  a  ‘linearizing  effect’ 
on  the  minimization  of  E(a,b)  without  losing  much  accuracy.  Keeping  in  view  the  above 
intuitive  idea,  the  modified  error  Era(a,b)  is  defined  as 
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where 


(3.11) 


Em  (a,  6) 


i=ifc=iy"** 

—  (g(joj)h(joj)a(ju))dlk{jut)  +  n»fcOu,))l  ^ 


-■  ^2  atrij^crkU^)  > 

r--l 

n 

rnlk{ju;)  =■-  hiT{ju)dTk{ju)  , 

r_  1 
n 

nlk(]*>)  -  /ir(ja>)mr*(;u;)  ,  l<i,fc<n 

r-  1 


(3.12) 


The  above  ‘linearization’  approach  for  SISO  system  design  and  model  reduction  has  been 
effectively  used  by  several  researchers,  including  Levy  [10],  Rao  and  Lamba  [20],  Rattan 
[18],  Rattan  and  Yeh  [19],  Pujara  and  Rattan  [15,16,17],  and  Pujara  [13,14].  For  the  sake 
of  brevity,  from  now  on,  we  will  write  Era(a,b)  as  just  Em. 

Before  we  begin  to  take  the  partial  derivatives  of  the  error  function  Era  with  respect 
to  the  unknown  controller  parameters  and  set  it  equal  to  zero,  we  need  to  simplify  Eni. 
First,  define  the  frequency  responses  of  the  (i,k)ik  entries  of  the  plant,  the  feedback 
element,  the  ‘desired’  and  the  controller  transfer  function  matrices  as  follows 


and 


where, 


glk{j *>)  ~  -4,*(u>)  ‘  u-)  ,  g(j^)  -  A{uj)  +  >#(->) 

htk(juj)  -  /■.',*( u;)  +  jFlk{ u--)  ,  Mju.’)  -  L(‘*'')  ■+  j^'(u') 

d,k(j^)  -  L,*(u;)  +  jA/.kfu;)  ,  <f(ju>)  L(u>)  +  jM{u) 


f(«» 

cxk(j<jj)  -  £  ai,,fc)(j^)r  •"  *.*(«*>)  +  jftkM 

r  =  0 

a(;u>)  =  1  +  5Z  Mi4*')”  ■-  <r(v)  +  i^(w) 

r-r  1 


(3.13) 


(3-14) 


(3.15) 
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=  4''k)  -  4''k)"3  +  a<?k)u>*  -  +  . . .  ) 

A*M  =  a^a,  -  a<**>u,3  +  -  4*»u,7  +  •  • . 

v(u>)  =  1  -  62u>j  +  btu> 4  -  66w®  +  •  •  • 

0(u>)  =  ijw  -  63u>3  ■+  6su>8  -  67w7  +  •  •  • 


Using  the  above  notation,  we  define  the  following: 

V>H  =  R{fc(jw)rf(iu»)}  =  E(u>)L{u>)  -  F(u>)M{u>) 
9(v)  =  S{h{j*)d{ju>)}  =  E(u>)M(u>)  +  F(u)L(  u>) 
<f>(u>)  =  R{£(ju>)Mju’)}  =  —  B(u>)F(ut) 


and  also 


7(u>)  =  3{5(ju>)/i(ju>)}  =  A(u)F(u>)  +  B(e>)E(u>)  j 


Q.*M  = 


*ik(u)  = 


*l{lik(jv)}  =  R  |^^r(j'u;)cr*(jiu;)| 

n 

£[4„(u>)crr*(a;)  -  Bir{ u>)0rk(v)} 

r=l 

S{fo(iw)}  =  9  |5^^r(j'u;)cr*(i7'u;)  j 

n 

=  +  ^.r(u>)<Tr*(u;)] 

Til 

=  ^(£tr(u,)Ir*(u,)  -  Fir{u>)Mrk(u,)} 

T  =  1 

=  3{m,*(;u>)}  =  3|£fc<r(>w)d;*0w)J 

=  fl\Eir(u>)MTk{»)  +  FfrHM*)] 


fal 


(3.16) 


(3.17) 


(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 
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where  *{•}  and  /{•}  stand  for  real  and  imaginary  part  of  a  function,  respectively.  Now 
finally  to  simplify  the  modified  Error  Era,  we  let 


-  %{:h{ju))d(ju>))llkUu>)} 

(3.26) 

SR{iv,('*0  4  J0(u;)][^,iit(u;)  +  >/*.*(“'))} 

(3.27) 

=-  *,'>(«*>  iQ.k(->)  - 

(3.28) 

=  V;(v(w).4tr(a>)  -  «(w)Bir(w))arfc(w) 

r-cl 

-(0(*M.r(w)  +  rp(u>)Bv  ,(u>))/Mw)l 

(3.29) 

-  ${lh(ju>)d(ju>)}llk(ju)} 

(3.30) 

3{  V’(-‘)  +  jtf(u’)l(Q.fc(u>)  4  J  } 

(3.31) 

=  $>,*)<?»*(“>)  -f  ^(u))A,fc(a-) 

(3.32) 

0(u;).4„.(u>)  4  il>{u)B,r(u}))(Trk(v) 

-  1 

-r(l/>(u>)/lir(u>)  -  ^(w)F,r(^))/?rfc(w): 

(3.33) 

=  ft  {  g !  j )h ( jiu>) j Q (.?u>)<f,fc ( ju ) } 

(3.34) 

=:  ft{  c(--)  -  j7(u)j|<T(w)  4  j(3{ui)}[Llk(^)  4  j  A/,  *(“■').} 

(3.35) 

-  ;(^(-?)£lfc(w)  —  -r(w)A/i*(u^))cr(i^) 

-  (d>(u>)A/,*(u>)  4  7(w)Li*(w))/3(w)] 

(3.36) 

=  0{|ff(jw)/i(ju;)iQ(jw)dlfc(;u»)} 

(3.37) 

=  ${^(«0  +  >7(^)]k(w)  +  30{  “>)][£.fc(w)  +  jA/,fc(u>)]} 

(3.38) 

=  [(<£{u/)M,fc(w)  4  7(u/)£,*(u/))(r(u/) 

4(^w)i,ii(w)  -  7(w)Mi*(u»))/3(u;)] 

(3.39) 
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A\*(u>)  =  R{nlk(ju>)}  -  ft 

=  *{£[(<?-(")  +;iitr(a,))(Vv^)  4 

= 

r-l 

-•  ]T  Y,  [(^ir'(w)V’rfc(w)  --  Bir'{ii>)0Tk{u))ar.r(u>) 

r  =  l  r'-l 

Vifc(w)  =  ^s{nik(ju>)}  -  3  |5Zilr(ju,-)mrfc(;u;)| 

-  4-  j/iir(u.-))ftv*(-0  + 

-  £[$„(*)*,*(*)  +  ^,(u.-)v,*(w); 

r-.l 

-  ^  ;(.4,r»(u>)0,*(t*.’)  T  5,r'(-^‘)Vrfc(^))0’,.r(4-’) 

Til  r'il 

-(.4,,.( Ul)ll>rk(uj)  -  B,r‘{u))8Tk(u)))0r'r{u) 


(3.40) 

(3.41) 


(3.42) 


(3.43) 


(3.44) 

(3.45) 


(3.46) 


(3.47) 


Substituting  Equations  (3.26),  (3.30),  (3.34),  (3.37),  (3.40)  and  (3.44)  into  Equation  (3.11), 


we  obtain 


E’n  =  it  it  f  4  KM.kH  +  j\k{u)} 
.=1  1 


{{Ttk(*)  4-  jVik(»))  -  (Xtk(u>)  4  jYlk( u-0)}i2  d* 


(3.48) 


=  E  E  /""  {(«.*(")  -  *»(*>)  -  r„(w))1 

*-.i 

4-(A,fc(u;)-y;jk(w)-t/.fc(u;))J}^ 


(3.49) 


Now  we  minimize  the  error  function  En  in  Equation  (3.49)  by  equating  to  zero  the  partial 
derivatives  of  Em  with  respect  to  a[^  and  bz,  for  0  <  r  <  p(i,k),  1  <  i,  k  <  n,  and  1  <  z  <  q. 
It  will  be  shown  that  this  will  lead  to  linear  algebraic  equations  involving  those  unknown 
controller  parameters.  For  the  sake  of  brevity,  from  now  on,  we  will  express  a  function  of 
ui  just  by  its  function  name  (i.e.,  f  =  f(w)).  Now  we  divide  the  partial  derivatives  of  Em 
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with  respect  to  those  unknown  controller  parameters  into  four  major  partial  derivatives  as 


follows: 


dEm  _  0 

8avw  “  ’ 


r  -  0,2,4, .  . .  ,p{v,u>)  , 


1  <  v,u'  <  n 


(3.50) 


5£m  df3vw  Q 

r  =  1,3,5,...  ,p(t>,tu)  , 

5aiVtW) 

aa/,u') 

1  <  r,  v>  <  n 

(3.51) 

5£m  da  _  _ 

2, 4, 6, . . . , 9 

56, 

5cr  8br 

(3.52) 

5£m 

56/ 

8Em  83 

83  8br  ’  r  ‘ 

1,3,5,. ..,9 

(3.53) 

In  order  to  do  the  partial  derivatives  of  Em  as  in  Equation  (3.50),  we  need  the 
following: 


where 


1 1  p  f2(w  -  Xu,  -  r.o  ( pk  -  p  -p] 

,r ]  k-l  ^  u,»*  [  \0(TVW  (y(Tvw  d(Tvw  J 

.  2f>  Y  „  dY,k  dV,k 

+  -  >,*  -  Uik)  a - -  7T - -  5 - 

y  v^vu.  Cr(TVUJ 


(3.54) 


g** 

5<r„u. 

-  V^i«  -  , 

—  Btv0wk  > 

f't/ VUI 

7T^  0 

oavw 

(3.55) 

-  +  vB,v  , 

-  A,v6wk  +  Btvipwk  > 

jLU*  -  0 

ocvw  t 

i  vill 

/  5/1, fc  5A\*\ 

(3.56) 

= 

y  ^O'vw  / 

(3.57) 

— 

t ft A\v  & $iv  A{v^Pwk  “f“  i?,A* 

/  5A,*  5K,*  \ 

(3.58) 

new  _ 

Bxk  ~ 

\0tfvu>  J 

— 

6Aiv  +  $Biv  -  Atv6wk  -  £** 

rtu .4. 

(3.59) 

So  that, 
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,  X  T  )( *£±  _  = 

±[MW^  -  6B")°rk  -  M?(eA"  +  1>Bir)P ,fc) 

r=l 

—  £  £  {A^iMrk  ~  Btt.6rk>r‘r  ~  A^(A,r.6rk  4  #,'VV*)/U] 

r=l  r'=l 

-  iMik)*  -  A^(4>Mik  4  nLik)fi] 


(3.60) 


(Ai*  -  Kf*  -  l/,*) 


aA,t  _  ev^N  = 


j^[B^(6Air  4  ***>,*.  4  B^Atr  -  6Btr)(3rk) 

r=l 

—  ]§P  +  ^tr'V’r*)<,V'*  4  (■'^tr'Vv*  ~  -S»r'0rfc)/?r'r] 

f  =  l  r’=l 

-[Br(Wfc  4  7Ilt)<r  +  -  7M.fc)/3] 


(3.61) 


Then 


(/*.*  -  AT,*  - 


ar.t\ 

VU>  J 


=  £[(lM?r  +  gBt7)Air  -  (tA?  -  rl>B^)Bir]<rrk 

r=l 

-  tl(M7  -  +  (V-47  +  liDH'IA* 

r=l 

-  £  Eitfsv* + *»*»)*-•  -  (irr«.*  - 

r=l  r'=l 

+  E  E  K>isr^fc  -  B^MAi  r.  +  (isrvv*  +  Br^fc)^.]/3r.r 

r=l  r'=l 

-  7  A/.*)  4  B;?{4>Mik  4  7^*)]* 

-[BS“(^*  -  7Af<0  -  iS“(^M<*  4  7li*))^ 


(3.62) 
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(3.63) 


where 


and 


=  £  -  £  *SA.  -  £  £  jr. +  f;  £  iS,.A- 


f=l 


r=l 


r=J  »<=1 


r=l  r'=l 


-  N"0 


fVU) 
ikr 


l 

Jvxv 
Jxkrr • 


IS  vu> 

A  ikr 


■■  CpA„-DXBir 
W£Air.-V£!Bi,. 
i>™Atr  +  C°?Btr 
likrr‘  =  V^Atr.  +  W,^Btrl 
A??  A*  +  B^Bik 
BifAik  -  A”?Btk 


Msr 


-  -  -rMik  , 

=  <f>Mik  +  7^ifc 

=  iMJj*  +  6B"?  , 

=  $A£-1>BS’ 

=  isrvv*  +  Bsr$rk , 

kst 

=  AXOrk-BW  rk 

Now,  substituting  Equation  (3.63)  into  Equation  (3.54),  we  obtain 


(3.64) 


(3.65) 


8E„ 

da^ 


££ jT“ 2  E  ft**-  -EE 

«=1  fc=l  •»  lr  =  l  r=l  r'=J 

-  £  Asia.  +  £  £  is,. a* 

,  =  1  r=lr'=J 

-Jtirr*  -  ivy*]  **  ' 


(3.66) 
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n  n  n  ,u>  n  n  n  n  fU' 

E  E  E  /  “  -  E  E  E  E  / 

-  E  E  E  f + E  E  E  E  f" 

iZ 1  fc=l  r=J  ■'*■'•*  i=U=lr  =  lf'=rw‘» 

-  E  E  r'“  -  E  E  /":‘ 

i:  1  fc=l  »=1  *=1  ** 

n  v.  (  n  r  .  w»  n  />wr  ,  . 

ee2Ie  /  ” -  e /„ •; 

-  E  [ /“:‘  -  E 

r  =  l  rw'»  r'=l 


(3.67) 


-  f  "  -  [/^‘  *7/^  } 


(3.68) 


n  n  n  -wf  n  n  n  #w'. 

EEE/'C'-^  -  EEE/' 

^k=U=l^‘  ,=H  =  1  i=l 

n  n  n  -w' 

=  EEE  /  ’*  &»»*■> 


»=1  *  =  17  =  1^* 


by  letting  r  =  i  and  i  —  r 


n  n  n  n  n  n  n  n  f  a/  .  „ 

E  E  E  E  /  -  E  E  E  E  / 

»<  =  lr  =  l  »=1  fc=l  “'•* 
n  n  n  n 

=  D  £  £  £  /  ~ 

i  =  l  *  =  1  r  =  l  r’=l  ■,w”' 


by  letting  r'  =  t,  r  =  fc,  t  =  r  and  k  =  r' 


Finally,  putting  it  all  together  as  in  Equation  (3.50),  which  is 


&  Em  d<Jvw 
d<rvw  da[v‘w) 


z  -  0, 2, 4, . . .  ,p(v,iu)  , 


1  <  v,w  <  n 


(3.69) 


Therefore,  we  obtain 
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UAi  1C  -  %X: 

- 1  [/:;  - 1 1  Sa**** 

-[C^^H-IjCj^H} 


z  =  0,2,4,.  ..,p(v,w)  , 


1  <  v,w  <  n 


(3.70) 


which  implies 


where 


n  n  (  n  f  ,w'  Q_  n  ,w’  , 

s  s  te  [/..;  -  s  4: 

-£  X:  &>****-&£ 

n  n  fl,r  3  =  °-2’4 . P{v,u>)  , 

=  EZ  /  ""ZsLMS’du  , 

J  <  U,TT  <  77 


(T  -  (7  1  —  — 62^^  ~f~  b\U)*  -  6gUJ^  -{-  •  •  • 


(3.71) 


(3.72) 


As  an  example,  take  z  =  0  and  fix  a  value  of  {v,u).  Note  that 


0(TVW  O  .  (v,u;)  (v,u;)  2  .  1  1 

- 7—  "i  =  - 7 - r  An  “  Oj  ^  rr  1 

o  (f.u*)  Q  (ir.tu)  1  u  *  J 

da[0  day  1 


After  some  algebraic  manipulation  Equation  (3.71)  leads  to 


(3.73) 
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We  now  use  the  following  set  of  notations 


(3.74) 


cO'.U') 


(3.75) 


Since  a ik,  /?ik,  o  and  /?  are  given  as  in  Equation  (3.16)  and  using  the  above  notations, 
Equation  (3.75),  we  can  rewrite  Equation  (3.74)  as  below: 

E  E  [To,’*’v,u',Go  ,fc)  -  T^Kv'w)a[l'k)  +  •  •  •] 

k=i 

-EE  -  Ui,Xv'w)a[''k)  +  •  •  •] 

«=i  *=i  1 

-  [-5^u'“’)62  +  S[v'w)b< - ] 

-  [fl(1v’w)61  -  R[v'w)b3  +  ■••]  =  S^lU,) 

(3.76) 

We  can  similarly  generate  equations  similar  to  Equation  (3.76)  by  taking  partial  deriva¬ 
tives  of  Era  with  respect  to  o^V,u>^  for  z  =  2,4,6. ..,p(t/,w)  and  1  <  v,  u>  <  n.  These  equations 
can  be  written  in  a  matrix  form, 

XiYi  -  Zi  (3.77) 
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where  the  matrix  Xj  is  given  by  Table  3.1 


rp(i,k,V,w)  y(i,k,v,w)  rp(  l,k,v,w) 


p(v,w)  _  ^(v.ui)  _ 


(3.78) 


Table  3.1:  The  matrix  Xi 


Y,  = 


(3.79) 


Now  we  establish  similar  relationships  as  above  by  doing  the  partial  derivatives  of 
Em  as  in  Equation  (3.51).  First,  note  that 


(3.80) 


&Hik 

dXik 

dTtk 

dPvw 

dpvw 

d@vw 

dxtk 

dYlk 

du%k 

,dpvul 

dpvw 

dpvw 
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where 


=  -(8Aiv  4  m.)  ,  +  B^wk)  ,  |£-0 


=  xpAlv  -  6Blv  ,  -  A,„Vw  -  0 


(3.81) 


Now, 


dmk  _  dXjk  \ 
#&„,  80vw ) 
8Xlk  _  8Ylk  \ 
80vw  80vw) 


-{0Aiv  +  rf>Biv  -  Aiv9wk  -  Btvipwk)  =  -By? 
lM,„  -  $Blv  -  Aivyl>wk  +  B,vBwk  =  A*? 


(3.82) 

(3.83) 


Thus,  we  have 


(*x.fc  Xik  T«)(d/3vw 


8Xik\ 

0/W 


f^l-B^Air  -  6B„)oTk  +  B:?(9Air  +  rl>B,T)0rk} 

r  =  l 

+  ]T  Y,\B??{Mrk  -  Bir.6rk)<TT'r  ~  B^{AXT,6Tk  +  Bir^rk)0r'r] 

r=l  r'=l 

+|B3T(*i.»  -  lM,k)c  -  BZ^Mik  + 


(3.84) 


(A„  Y,k  U,k)  (^‘  s^)- 

E(4S-(M„  +  V'B.-K.  +  AXWA,,  -  } 

r=l 

—  [AJjJl’(Air>^r*  +  ^r<V,ffc)o’r'r  +  AJjJW(  Air'Vrk  ~  B,T<6 rk)0TiT] 


r=l  r'=l 

-1  A^{d>Mik  +  iLik)o  +  -  7M,fc)/3] 


(3.85) 


Then 
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fe  -  S + -  *w  (sc "  © 

=  D(«zi~  -  Mr)*.  +  (Mr  +  «B.r)B..]<’.» 

r=l 

+  EKimst  +  «*srM.r  -  (Mr  - 

r- 1 
n  n 

-EE  -  B™ri>rk)Airi  +  Urrvv*  + 

r=l  r'=l 

-EE  lUaVv*  +  b,7M^>  -  (irr^*  - 


r=l  r'=l 

+\W?(<t>Lik  -  1Mxk)  -  A*n 4>Mxk  +  7M]* 
-  7^<*)  4  4  7 Lik))0 


(3.86) 


E  *.*>'*  +  E  -EE  ^SKr-^v  -EE 

r=1  r=J,<=i  r=lr'=J 

+*3r<r  -  M.7/3 


(3.87) 


by  using  the  same  notation  as  in  Equations  (3.64)  and  (3.65). 

As  before,  by  substituting  Equation  (3.87)  into  Equation  (3.80)  and  by  using  the  argument 
similar  to  the  one  used  earlier,  we  obtain 

fir  n  n  f  n  f  n  MUf> 

|f=  =  El!l/"«^-Erc.^ 

Opvw  <ss n_i  U=1  L-'w-*  r«= 

r=l  rw-»  r'=l  J 

+  -  [£>~H} 

Finally,  putting  it  all  together  as  in  Equation  (3.51),  which  is 

BEm  _  8Em  60,w  *«l,3,5,...,jp(t»,w), 

dQ(1v'“)  0a<*.»>  ‘ 0  ’ 

1  <  v,w  <  n 


(3.88) 


(3.89) 
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Therefore,  we  obtain 


obtain 

£  s 2  {£  [C  - 1,  £r 

+£  \C  - 1,£: 

+[£‘^H-  [£‘£feA>Hl 


2  =  1,3,5, .. .  ,p(v,w) 


which  implies 


t',U'  <  71 


.  . *  J 

.  1  f  t!,w<n 

where  a  is  given  by  Equation  (3.72). 

As  an  example,  take  z  =  1  and  fix  a  value  of  {v,u).  Note  that 

a 


da[v'w)  da!; 


After  some  algebraic  manipulation  Equation  (3.91)  leads  to 

fc=l  Ir^l  r'=  \Jwrr’  / 

+£[(/*-*=- 

rrllV“-‘  r'=l / 

+  [(£‘“^H-K£‘^H} 

«  n  »w'  1 

-EE  / 

i=U=l  J 


(3.90) 


Ss  is  1C  $***«•*■ - 1  c 
-  £  tC  £='•■“-  - 

-f  <^t,u'  jirvu-  •  i  5,3,. .  ) 

f  k  &?=»*“  H "  IL  5jS«w*-j} 

a_ff[  /“I.  «A.  -•  -  ••*.* . **,*) .  (*•«) 

£ri&V~*  “  j  ’ 

1  <  V .  W  <  71 


(3.92) 


(3.93) 
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By  using  the  same  notation  as  in  Equations  (3.75)  and  using  ait,  /?ik,  o  and  /?  as  given  by 
Equation  (3.16),  we  can  rewrite  Equation  (3  93)  as  follows 


£  £  [v «  +  ...] 

izl  k=l 


i=i  k=i 


4  \-R[v  w)b2+  R{v'w)b<  -  ] 

-  [Sf^,  -  +  ■••]  =  - M-1 


(3.94) 


We  can  similarly  generate  equations  similar  to  Equation  (3.94)  by  taking  partial  deriva¬ 
tives  of  Em  with  respect  to  c^’^for  z  =  3,5,7,...p(t/,4/)  and  1  <  v,  u  <  n.  These  equations 

z 

can  be  written  in  matrix  form, 

XjYj  =  Zt  (3.95) 


where  the  matrix  X2  is  given  by  Table  3.2. 


rp[i,k,v,  u) 


ui'Xv'w) 


-S<UlU)  -R[v'u)  s<ViU,) 


—  u> 


rp(i,k,v,\v) 


£/(*  ,k,V,w) 


s{v,u’) 


(3.96) 


Table  3.2:  The  matrix  X2 


and 
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,  z2  = 


(3.97) 


r[v-w) 


Now  we  establish  similar  relationships  as  above  by  doing  the  partial  derivatives  of 
Em  with  respect  to  b’s  as  in  Equation  (3.52).  First,  note  that 


where 


V'  V  T;‘  f 9f  -  x  T  \  ( dti,k  dXlk  -  — 

,  ,  (d\,k  d\\k  dVikW  . 

+  2(A,fc  -  U  V,k)  -  5ct  "  da)\ 


(3.98) 


_  0  ]  .  o  ^  ~  (f,Llk  -  lM-.k  =  'U 

_  n  Q£xk.  _•  0  , 

'  1  ’  ITo-  '  OCT 


Thus,  we  have 


(3.99) 
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(3.100) 


-  j2[AM’A>r  -  0B„)<rrk  -  Aik(0Ai,  4  tl>BtT)(3Tk] 

r=l 

4  ^2  53  [AikiAir'Iprk  ~~  Bxr’Qrk^r'r  ~  Aik(AtT'8Tk  4  Bir'  )@r'r] 

r= 1  r'—l 

4 [Aik(4>Lik  -  'iMlk)o  -  Aik[<t>Mik  +  7 Uk)P) 


—  "I"  ^Bir^ffrk  4  BikfyAir  ~  &  Bir)#rk\ 

r=l 

4  53  53  [^»fc(^»r*^rfc  +  Bir'lprk)0*' r  +  Bik{AiT>yprk  ~  Bir'0rk)Pr'r\ 

*■=1  r'=l 

4[5ifc(^Mi*  4  +  *ih(<f>Lik  ~  *yAfi*)^] 


(3.101) 


Then 


-  **  -  r“>  (-^) +  -  *  -  «w  (-^) 

n 

=  ~  £1(^4*  +  6B,k)A,r  -  (£4,*  -  V’5,*)B,r]<Tr4 

r=l 

n 

+  £l(**fc  ~  ypBik)Air  4  {ipAik  4  0£i*)£„]/?r* 

r=l 
n  n 

"l"  5  ^  5  {{Aiklprk  4  BikQrk)Air>  (^4|*^r*  *ik^Prk  )Bir>](TriT 

r=l  r'=l 
n  n 

—  53  53  !(***'*  “  Bik^>rk)Air<  4  (Aiklprk  4  I?ifc0rfc)i?ir>]/?r>f. 

*=1  r'=l 

4[-4,*4,*  4  BikBik)<r  4  [B,*y4i*  -  A***}# 


(3.102) 


n  n  »  "  _ 

=  -  53  likrOrk  4  53  ^ikrPrk  +  51  53  ^ikrr'^r'T  ~  53  53  Likrr'Pr'r 


»=] 


r  =  l 


r  =  J  r'=l 


r=  1  »'=1 


4Af,*<7  4  A\fc/0 


(3.103) 
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where 


hkr  —  CtkAxr  -  D,kBxr 
Jxkrr’  -  VVxkrAlT.  -  V,krB,r. 
■ft'.fcr  =r  DlkAxr  -f  CxkBir 
Likrri  —  VikrAiT<  +  WtkrBxr< 
Mxk  =  AxkAxk^BxkBxk 
A’.*  =  BxkAxk  -  AlkBxk  =.  u 


(3.104) 


and 


Aik 

-  <t>Lxk  ~  lMxk  , 

Bxk 

=  4>Mxk  -  •'I,* 

cik 

=  i’Aik  +  6Bik  , 

Dxk 

=  &Alk  -  vBlk 

> 

(3.105) 

wxkr 

=  Aiklprk  +  Bik9rk  , 

^.kr 

~  Alk9rk  —  Bxtvrk 

As  before,  by  substituting  Equation  (3.103)  into  Equation  (3.98)  and  by  using  the  argu¬ 
ment  similar  to  the  one  used  earlier,  we  obtain 


JL  f  "  f  fu‘rk  n  : 

d(T  E  £  { - £  ^  irkx<TXkdu  -  £  j„.klaxk<L>^ 

£  \j  K*kx0xkdu>  -  ^  J  Zrr'k, 0ikdw 

[£‘ + 

Finally,  putting  it  all  together  as  in  Equation  (3.52),  and  noting  that 


+ 


+ 


(3.106) 


dE^ 

dbt~ 


dEm  da 
da  dbt 


=  0  ,  z  =  2,4,6, ... ,q 


(3.107) 


We  get, 
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££*{-£  \r  -  i  /"”■ 

t  —  1  /t~- 1  l  r=l  l/"*»  °'6*  r.=  1-/u.rr,  yOj 

V-'  17  ,o  J  V"  /“'«■'  -y- 

+  fr:‘  1 

•'**>,4  C/O* 


/  "  —  A \kpdw 
•-.»  dbr 


=  0,:  -  2.4.6,. ..,g 


which  implies 


n  n  (■  n 


ssi-siri7-— -sci7— 

di.K'kAk,L,~%JZ 

,  5vw-“H  -  [/...  3^‘H } 

a  "  *r  /-■:»  a<r__  i 

■“"■2-2-1/  ,  z  =  2,4,6,...,9 

»-i  k=i  !/-•*  O0*  J 


C'  0CT 


+  E 

r=l 


L 


where  a  is  given  as  in  Equation  (3.72). 
As  an  example,  take  z  =  2.  Since 


It  can  be  easily  seen  that  Equation  (3.109)  leads  to 


=  —U) 


EE{E  (/  ^Irki-  J2  J  r  v2Jrr,kl\  aikdw 

i  =  l  *  =  1  V.  r—  1  V"*»  r,_,  /w,,,  / 

{[y‘K.„-±['"^L . )a^ 


-E 

r  =  l 


'*  u>2 <7<k>  -  '*  u/aJV,k)  j3d.j 


y  u)3Mtkdw 

We  now  use  the  following  set  of  notations: 


-EE 

i-  1  *  1 


(3.108) 


(3.109) 


(3.110) 


(3.111) 
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(3.112) 


Using  aik,  Pn,  a  and  P  as  given  by  Equation  (3.16),  we  can  rewrite  Equation  (3.111)  as 
follows: 


«~1  k=l 
n  n 


-EE 

«=i  fc-i 

—  [— Z462  4-  —  •••]  =  Zj 

(3.113) 

We  can  similarly  generate  equations  similar  to  Equation  (3.113)  by  taking  partial  deriva¬ 
tives  of  Em  with  respect  to  bz  for  z  =  4,6,8,...,q.  These  equations  can  be  written  in  a 
matrix  form, 

XSY,  =  Z, 

(3.114) 

where 


X3 


P?'k)  -Q{x'k) 


-pr]  Qi 


-P4(a)  q['a) 


^4 


pW 


(•» 


0  -Zg  0 


(3.115) 


Table  3.3:  The  Matrix  X3. 


30 


,  Z3 


(3.116) 


Now  we  establish  similar  relationships  as  above  by  doing  the  partial  derivatives  of 
Em  as  in  Equation  (3.53).  First, 


where 


dfiik 

dXlk 

dTlk > 

,  d0 

dP 

0/Ty 

(d\ik 

dVik 

V  dp 

dP 

*0 

(3.117) 


=  ~(Wi  +  -ri,»)  -  -B* 


vft  '  ~df  ^ 0  ’  - 7.w.fc  - /U 


Thus,  we  have 


(3.118) 
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(3.119) 


(,,-A fc-7M(-f£)  = 

n 

-  9Bir)crrk  -  Bxk{9Axr  +  i>B„)3rk ] 

r  =  l 

n  n 

~  XI  XI  -  B,r,0rk)crrlr  -  Bik(Air.0Tk  4-  i?,,<Vrfc)3r'rj 

r  =  l  r'=l 

-[J3,-* (?>£,*  --  'yMtk)cr  -  Bik(4>Mxk  +  7i,*)/3j 

-  £{2, •*(**„  +  V-#.*  )<*••*  -f  AiktyAir  --  8Bir)9rk] 

r  =  l 
n  n 

+  XI  XI  Bir<lfiTk)aT>r  +  Alk(AiT'^rk  —  BlT>0Tk)flr'r) 

r  =  l  r'  =  1 

+  [>M</>M,*  +  iLik)cr  +  Alk{<i>Lik  -yMik)0) 


Then 

(**  -  **  -  Tik)  )  +  (A*  -  V,fc  -  IT*)  (-^) 

=  -  £[(«*  -  +  WAik  +  9Bxk)BXT]cTtk 

r=l 

—  XlKV’Afc  +  8Bxk)AXT  —  (0i4jk  —  rl>Bik)Bir}{3rk 

r= 1 

-f  XI  X/^'^‘fc^rfc  ~  Bikll>Tk)Air'  +  {Aikfprk  +  Bxk9Tk)Bir']crriT 

r  =  X  r'  =  1 

+  y,  XI  [(•^‘fcV’r*  +  Bik0rk)Air'  ~  ( Axk8Tk  ~  Bxklf>rk)  #ir’]/?r'r 
r=l  r'=l 

-[BikAxk  -  AikBxk]cr  +  [AikAik  +  BxkBx k]P 


XI  K*r<Jrk  -  XI  IikrPrk  +  XI  XI  Likrr'(Tr,r  +  XI  E 

r;;;1  r  =  1  f=W'=l  r  =  lr-  =  l 

-Nikcr  +  Mik/3 

by  using  the  same  notation  as  in  Equations  (3.104)  and  (3.105). 


As  before,  by  substituting  Equation  (3.122)  into  Equation  (3.117)  and 
argument  similar  to  the  one  used  earlier,  we  obtain 


(3.120) 


(3.121) 


(3.122) 


using  the 


32 


~  -  £  £ 2  {  -  £  f  £  LMkiu 

0>D  t=l  *-  I  l  r=l  r"'*  r,._zl  J^rr' 

-  E  \r 4  -  E  Jrr'k,3lk(Lv 

Finally,  putting  it  all  together  as  in  Equation  (3.53),  and  noting  that 


^  ^  0  ,  z  =  1,3, 5,...,  9 

dbz  d(3  3bz 


Therefore,  get, 


e  e  2  ( -  e  r-  —-Kr^du  -  £  r'rri  ^-Lrr.kl<TXk<L, 

ht,  l  db>  ,tV“W  db* 

"  r:»  a/?  "  t<r.  d(3 

-E  /  AirlrklP'kdU;-  E  /  ^T-Jrr'ki&kdu 

r- 1  l/wrfc  r»-l  7wff'  vbt 

-  £  + [£‘  } 


=  o  ,  z  =  1,3,5,. ..,9 


which  implies 


§  S  {  5  [■£  at!  -  £  £" 

-EE  =o.*  =  1.3, 5...  .,9 

k=l  IA-‘  db‘ 


where  <7  is  given  as  in  Equation  (3.72). 


As  an  example,  take  z  =  1.  Since 


(3.123) 


(3.124) 


(3.125) 


(3.126) 


aa  a 


at,  96, 


-  [t,uJ  —  630?^  ~r  ‘  ‘  ■  ]  -r  uJ 


(3.127) 
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After  some  algebraic  manipulations,  Equation  (3.126)  leads  to 


«=1*  =  1  l  r:l  1  x-1--*  r,_j  ) 

'  5  ((/.. uU 

n  n  r  /u/r 

=  E  £  /  “  -  o 

i=i  *=i  l/".* 

(3.128) 

By  using  the  same  notation  as  in  Equation  (3.112)  and  using  /?ik,  a  and  /?  as  given  by 
Equation  (3.16),  we  can  rewrite  Equation  (3.128)  as  follows: 


-EE[«SwMa'-<?W  +  -] 

«=ifc=i  J 

-EENaMa,-«ia,4“'  +  -l 

+  [Zjij  -  Ztb3  +  •  •  -J  =  0  (3.129) 

We  can  similarly  generate  equations  similar  to  Equation  (3.129)  by  taking  partial  deriva¬ 
tives  of  Em  with  respect  to  bz  for  z  =  3,5,7,.. .,q.  These  equations  can  be  written  in  a 

matrix  form,  where  (3.130) 

X4Y4  =  Z4 


q(.»  qO.*)  p0'fc) 


Z2  0  -Z4 


Q(>,k)  />(•.*) 


~Z4  0  Zb 


(3.131) 


Table  3.4:  The  Matrix  X<. 


V4 


a 


(•>> 

o 


c 


<*» 

2 


a 


'->) 


61 


0 


62 


63 


(3.132) 


Combining  all  the  matrices  in  Equations  (3.77),  (3.95),  (3.114)  and  (3.130),  we  can 


regroup  them  in  the  following  matrix  form, 

XY  =  Z  or  Y  =  X*Z 
where  the  matrix  X  is  given  by  Table  3.5. 


rp(l,k,V,w) 

_y(i,k,V,w) 

,jy(i,k,v,w) 

s{v'w) 

rp(i,k,V,W ) 

1  2 

_u(i,k,v,w)  .. 

—  j?3l’’U) 

c("'u) 

°4 

~  1  2 

J,(t, k,v,w) 

•  R[v'w) 

.  Ji.k.v.w) 
~  1  3 

rp(itk,VtH>) 

~~  1  4 

y(i,k,V,w) 

C(v.w) 

d4 

-r> 

1 

V 

•  • 

•  ^2 

0 

-Z4 

r^k) 

_pb>)  .. 

0 

0 

<?<•*> 

pO.M 

-4?(5>)  •• 

•  -^4 

0 

Table  3.5:  The  matrix  X 


(3.133) 


(3.134) 
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and 


Y  = 


(*.«0 

ao 


bi 

bj 

b3 


,  z 


$6 


(t’.u) 


~*i 


(v,w) 


•  3 


0 

z3 

0 


(3.135) 


Thus  the  matrix  Equation  (3.133)  determines  the  unknown  controller  parameters  aV,,k;  and 
bz  for  0  <  r  <  p(i,k),  1  <  z  <  q  and  1  <  i,k  <  n. 

As  an  illustration,  for  a  system  of  order  n  =  2,  a  second  order  controller  with 
unknown  parameters  a’s  and  b’s,  is  given  as 


C(-)  = 


a*1'"  +  +  «lMV  4,’!|  +  «i'  ”»  »  4’-’V 

+  a<”>S  +  <4”’  -  «iMI»  +  «iWV 


1  -*■  b\3  +  bis 2 


(3.136) 


For  this  case,  the  matrices  Y  and  Z  are  given  by  Equation  (3.137)  and  the  matrix  X  is 
given  by  (3.138)  (Table  3.6). 
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3.4  Remarks 

Before  we  conclude  this  chapter,  we  wish  to  make  some  explanatory  comments  on 
the  technique. 

1.  For  using  the  technique  proposed  in  this  chapter,  it  is  assumed  that  the  given 
plant  is,  at  least,  marginally  stable  or  if  unstable,  it  has  already  been, 
somehow,  made  at  least  marginally  stable. 

2.  The  ‘desired’  closed  loop  transfer  function  matrix  has  to  be  mathematically 
formulated  from  the  design  specifications. 

3.  Observe  that  the  integrand  in  Equation  (3.11)  for  each  fixed  i  and  k,  is  the 
difference  between  two  polynomials  in  ui  with  coefficients  involving  a’s  and 
b’s.  Note  that  when  we  equate  to  zero  the  partial  derivatives  of  Em  with 
respect  to  a’s  and  b’s.  we  will  get  linear  algebraic  equations  in  a’s  and  b’s. 

4.  Since  the  integrand  in  Em,  for  each  fixed  i  and  k,  is  the  difference  between 
two  polynomials  in  ui,  intuitively,  it  seems  logical  to  conclude  that  the  best 
‘fit’  will  be  achieved  when  the  polynomials  in  u  are  of  the  same  degree  unless 
the  higher  degree  terms  in  one  of  them  compared  to  the  other  have  relatively 
small  coefficients.  We  formalize  this  intuitive  idea  in  the  proposed 
algorithm. 

5.  In  order  to  check  whether  the  integrand  in  Em  for  each  fixed  i  and  k  has  the 
best  ‘fit’  as  described  previously,  we  can  just  check  the  pole-zero  deficiency 
between  each  entry  of  the  compensated  transfer  function  matrix  and  the 
‘desired’  transfer  function  matrix. 

6.  For  a  system  of  order,  n  =  2,  suppose  the  degree  of  the  plant,  controller  and 
the  ‘desired’  transfer  function  matrices  are  given  as  follows: 

8  {gik(s)}  =  mik,  1  <  i,  k  <  2 

Hg(s)}  =  n 
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where  mik  <  n,  for  1  <  i,  k  <  2. 

f>  {Cik(s)}  =  pik,  1  <  i,  k  <  2 
5{a(s)}  =  q 

where  pik  <  q,  for  1  <  i,  k  <  2. 

6  {dik(s)}  =  /ik,  1  <  i,  k  <  2 

6  {d(S)}  =  k 

where  ^ik  <  k,  for  1  <  i,  k  <  2. 

In  order  for  the  integrand  in  Em  to  have  the  best  ‘fit’,  one  would  like  the 
following  conditions  to  be  satisfied  as  far  as  possible: 

max  {mu  +  Pn,  m12  +  p2J  +  k 

=  max  {n  +  q  +  In, mu  4-  pi2  +  £21,  nii2  +  p22  +  £21} 
max  {mu  +  pi2,  m12  +  p22}  +  k 

=  max{n  +  q  +  l\i,  mu  +  pj2  +  £22)  m12  +  p22  +  ^22} 
max  {m21  +  pn,  m22  +  p2]}  +  k 

=  max  {n-f  q+  £21,  m2i  +  Pn  +  £u,  m22  +  p2i  +  ^11} 
max  {m21  +  pi2)  m22  +  p22}  +  k 

=  max  {n  +  q  +  £22,  m2i  4-  pn  +  £a,  hi22  +  P21  +  ^12} 

We  can  enhance  the  ‘desired’  transfer  function  matrix  D(s)  so  that  both  the 
polynomials  in  the  integrand  of  Em,  for  each  fixed  i  and  k,  are  of  the  same 
degree  by  adding  as  many  relatively  nondominant  poles  in  the  entries  of  D(s) 
as  needed  in  such  a  way  that  the  steady  state  response  stays  unaffected. 
This  approach  of  enhancing  the  desired  transfer  function  will  lead  to  better 
minimizing  of  Em  without  adding  to  the  order  of  the  controller. 
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8. 


If  the  coefficients  in  the  plant  and/or  ‘desired’  transfer  function  matrix  are 
large,  then  the  minimization  of  Era  will  lead  to  overflow  and/or  numerical 
errors  in  the  integration,  matrix  inversion  etc.  Thus,  to  avoid  this,  it’s 
advisable  to  divide  all  the  coefficients  of  the  transfer  function  matrix  by  a 
suitable  factor  to  make  them  small  for  computation.  Evidently,  this  has  no 
impact  on  the  algorithm. 

/ 

9.  The  frequency  interval  of  interest  [wik,Wik]  is  usually  the  bandwidth  of  the 
(i,k)ik  element  in  the  ‘desired’  transfer  function  matrix  D(s)  and,  in  general, 
one  gets  good  results  by  minimizing  Em  over  this  frequency  interval.  But,  in 
some  cases,  if  the  bandwidth  is  large  and/or  the  coefficients  in  the  transfer 
function  matrix  involved  are  also  large,  then  the  minimization  of  Era  over  the 
entire  bandwidth  leads  to  large  entries  in  matrix  X  and  vector  Z,  (Equations 
(3.138)  and  (3.137)),  respectively.  This  leads  to  overflow  and  numerical 
errors  in  the  integration,  matrix  inversion,  etc.  In  order  to  avoid  this,  it  is 
advisable  to  minimize  Em  over  smaller  intervals  than  the  bandwidths,  usually 
one  half  to  one  quarter  time  the  bandwidths.  One  may  have  to  tinker  a  bit 
with  the  limits  of  integration  for  Ero  within  the  above  guidelines  until 
satisfactory  results  are  obtained.  This  should  not  be  very  time  consuming  as 
the  algorithm  is  computer-aided. 

10.  A  computer  code  in  Fortran  77  has  been  developed  to  implement  the  algor¬ 
ithm  with  the  help  of  IMSL  routines.  The  program  has  been  written  to  be 
user-friendly.  The  listing  of  the  program  is  given  in  Appendix  A. 

11.  The  technique  is  very  easy  to  use  overall.  All  that  the  user  has  to  do  is  to 
input  the  parameters  of  the  plant,  feedback  element  and  ‘desired’  transfer 
function  matrix,  the  order  of  the  controller  (one  usually  starts  with  the  zero 
order),  and  the  limits  of  the  integration  for  Em.  As  output,  he  gets  the 
parameters  of  the  controller  and  also,  if  the  system  has  two  inputs  and  two 
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outputs,  he  also  gets  the  compensated  closed  loop  transfer  function  with  its 
poles  and  zeros  of  the  elements  of  the  transfer  function  matrix,  and  the 
frequency  responses  of  the  ‘desired’  and  the  compensated  systems.  If  the  ‘fit’ 
is  good,  one  is  done.  Otherwise,  the  user  increases  the  order  of  the  controller 
and/or  alters  the  limits  of  integration  for  Era  as  indicated  above  and 
continues  until  satisfactory  results  are  obtained. 

12.  Finally,  it  should  be  emphasized  that  this  simple  algorithm  would  not  solve 
all  the  problems  of  multivariable  control  systems  design  satisfactorily.  This 
will  only  give  a  candidate  for  a  preliminary  design  which  will  have  to  be  fine- 
tuned  to  make  it  workable.  But  it  gives  an  excellent  starting  point  for  design 
purposes. 

3.5  Summary 

In  this  chapter,  a  computer-aided  technique  of  designing  a  multivariable  control 
system  is  proposed.  The  controller  parameters  are  obtained  by  matching  the  frequency 
responses  of  the  compensated  closed  loop  system  and  a  ‘desired’  closed  loop  system  over  a 
certain  frequency  interval  of  interest.  The  technique  is  easy  to  use  as  the  controller 
parameters  turn  out  to  be  the  solutions  of  linear  algebraic  equations.  Finally,  some 
explanatory  remarks  were  given  of  the  algorithm  and  the  overall  technique. 
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4.  ILLUSTRATIVE  EXAMPLES 


Two  numerical  examples  are  presented.  First,  we  reexamine  the  example  of  Chen 
[1]  and  compare  the  controllers  designed  by  the  proposed  technique  and  Chen’s  technique. 
The  second  one  pertains  to  the  design  of  the  control  system  for  the  YF— 16  CCV  in  its 
longitudinal  axis. 


4.1  Example  1  [1] 

For  Chen’s  example  [1],  the  block  diagram  for  a  two— input  two— output  control 
system  is  given  by  Figure  3.1  in  Chapter  3  with  the  feedback  element,  H(s)  =  I2,  and  the 
transfer  function  matrix  of  the  plant  given  as 


G(-) 


1.265a  +  2.122  941.479a  -  1506.997 

7.190a  +  10.238  1241.0169a  -  2107.295 
a2  +  2.918a  +2.107 


(4.1) 


It  is  desired  to  synthesize  a  controller,  C(s)  so  that  the  overall  closed  loop  system 
matches  the  ‘desired’  closed  loop  system.  The  ‘desired’  closed  loop  transfer  function  as 
synthesized  by  Chen  [1]  to  meet  the  specification  is  given  by 


From  the  magnitude  of  the  frequency  response  plots  of  the  elements  of  D(s),  the  frequency 

/ 

interval  [u/ik.Wik]  =  [50,100]  for  1  <  i,  k  <  2  is  selected  over  which  the  error  function  Em  will 
be  minimized.  Using  Equation  (3.129),  the  zero  order  controller  is  obtained  as 
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Co(-)  = 


-212.5080  93.2877 


1.0392  -0.1393 


Therefore,  the  compensated  closed  loop  transfer  function  matrix  is  obtained  as 

F0  (»)  =  [la  +  G(j)C0(a)H(a)]-1G(a)Co(j)  (4-4) 

fn(»)  &(•) 

ft \(s)  fiiia) 

(4.5) 

where 

ftM  ~  709.5624(j  +  1.3102)(«  +  1.6084)(*  +  495.0238) 

/?,(-)  =  — 13.1391(j  -I-  0.9109)(a  +  1.3ll8)(«  +  1.6062) 

/ai(«)  =  — 238.268(a  -  0.0598)(#  +  1.3U8)(«  +  1.6062) 

/«  («)  =  497.865(*  +  1.3102)(«  +  1.6084)(«  +  704.6024) 
f°{*)  ~  («  +  1.3102)(«  +  1.6084)(«  +  485.5063)(*  +  724.8383) 

The  frequency  responses  of  each  element  of  the  compensated  and  the  ‘desired’  closed  loop 
transfer  function  matrices  are  given  in  Figures  4.1  through  4.8.  In  addition,  the  time 
responses  to  the  inputs  of  the  compensated  and  the  ‘desired’  systems  are  given  in 
Figures  4.9  and  4.10. 

The  zero  order  controller  designed  by  Chen  [1]  is  given  by 


K(3)  = 


-202.0704  96.3172 


0.9815  -0.1356 


The  compensated  closed  loop  transfer  function  matrix  is  given  by 
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Ffc(s)  =  [IJ  +  G(a)K(a)H(3r1G(3)K(3) 

fni3)  fti*) 

/?»(')  /»a(«) 

177 

where  fn(s)  =  668.4426(s  +  1.3102)(s  +  1.6084)(a  +  523.7637) 

/i2(-s)  =  -5.8233(3  -  0.006235)(3  +  1.31 18)(a  +  1.6062) 

/21(s)  =  -234.828(3  +  0.002073)(j  +  1.3118)(a  +  1.6062) 

/22(3)  =  524.2  3  8  8(3  +  1.3  1  02)(s  -+  1.6084)(s  -f  667.1694) 
f  (■*)  ~  (3  +  1*3102)(3  +  1.6084)(3  -f  516.830l)(j  +  678.7686) 

For  the  purpose  of  comparison,  the  frequency  response  of  each  element  of  Chen’s  compen¬ 
sated  closed  loop  transfer  function  matrix  are  superimposed  on  Figures  4.1  through  4.8.  In 
addition,  the  time  responses  to  the  inputs  of  Chen's  compensated  system  are  also  given 
in  Figures  4.9  and  4.10.  It  can  be  observed  from  Figure  4.3,  that  the  magnitude  response  of 
the  (2,1)  element  and  from  Figure  4.6,  that  the  phase  response  of  the  (1,2)  element,  of  the 
compensated  closed  loop  transfer  function  matrix  have  poor  match  with  the  magnitude  and 
the  phase  responses  for  the  corresponding  elements  in  the  desired  transfer  function  matrix. 
Therefore,  it  is  desirable  to  improve  this  particular  mismatch. 

Now  to  improve  the  match  with  the  desired  transfer  function  matrix,  we  consider  a 
first  order  controller  matrix.  Again,  using  Equation  (3.129)  with  the  same  frequency 
interval  of  interest,  the  first  order  controller  is  obtained  as 

I 

-224.879(3  -f  1.641)  99.618(3  +  1.575) 

1.00343(3  +  1.824)  -0.1212(3  +  1.831) 

Ci(3)  -  - - - - : _ ' 

(j  +  1.836)  (4.9) 


(4.7) 


(4.8) 
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Therefore,  the  compensated  closed  loop  transfer  function  matrix  is 


where 


Fi(a)  =  [I2  +  G(j)C1(5)H(S)]-1G(J)C1(3) 

fill3)  fu{3) 

fl  ii3)  flii3)  (4-11) 

7m>] 

/n(a)  =  660.236(3  4-  1.3102)(j  +  1.5494)(s  +  1.6084)(s  +  1.8223)(a  +  574.4794) 
f},(s)  -  -11.909(3  -  0.07757)(3  +  1.3118)(*  +  1.6062)(3  +  1.6131)(j  +  1.836) 

/ji (s )  =  -371.606(3  -  0.1277)(s  +  1.3118)(3  +  1.6062)(3  +  1.6608)(3  +  1.836) 

ft2(s)  =  565.8422(3  +  1.3102)(a  +  1.5494)(3  +  1.6084)(a  +  1.8222)(a  +  669.379) 
f'{s)  =  (a  +  1.3102)(j  +  1.5494)(j  +  1.6084)(a  +  1.8222) 
x(3J  +  1229.29623  +  380042.6328) 


Figures  4.11  through  4.20  show  the  comparison  of  the  frequency  responses  and  the  unit  step 
time  responses  of  the  compensated  system  and  the  ‘desired’  system.  These  figures  show 
that  all  the  magnitude,  phase  and  time  responses  of  the  compensated  closed  loop  system 
match  very  well  with  the  ‘desired’  closed  loop  system.  Comparison  of  Figures  4.1  through 
4.10  and  Figures  11  through  20  shows  that  the  results  obtained  by  the  first  order 
compensation  by  the  proposed  technique  are  much  better. 


4.2  Example  2 

A  preliminary  longitudinal  control  systems  design  of  YF— 16  CCV  airframe  is 
presented  using  the  frequency  matching  technique.  The  technique  is  used  to  choose  the 
parameters  of  the  controller  to  match  closely  a  ‘desired’  closed  loop  frequency  response. 
The  ‘desired’  closed  loop  frequency  response  is  based  on  meeting  the  military  flying 
qualities  specification,  MIL-F-8785C  [12,21],  for  Level  1  flying  qualities. 
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4.2.1  Airframe  Model 

The  airframe  model  chosen  is  the  YF-16  CCV  aircraft  in  conventional  flight  at  0.8 
Mach  at  20,000  feet  [7].  The  following  state  space  model  was  provided  by  Mr.  Tom 
Gentry,  the  Project  Engineer  who  derived  it  from  [7].  Based  on  these  aircraft  longitudinal 
data,  the  state  space  model  is  given  as  (Figure  4.21): 

x  =  Aiz  +  Biu  (4.12) 

n  —  rif  t  4-  n. 


where 


y  —  C  i  z  +  Dim 


A,  = 


-0.016484  111.35  -32.17  0 

-9.3359E-5  -1.0046  0  0.9983 


2.6976E-5  11.814  0  -1.1174 


6.648  8.662 


0.15306  0.16445 


21.096  4.9075 


0  10  0 


0  0  0  1 


(4.14) 


(4.15) 


(4.16) 


(4.17) 
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Figure  421:  Forces  and  moments  on  an  airplane  in  a  steady  climb 
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X 


u 

a 

e 

Q 


forward  speed,  ft/sec 
angle  of  attack,  rad 
pitch  attitude,  rad 
pitch  rate,  rad/ sec 


u  = 


S, 


horizontal  tail  deflection,  rad 
flap  deflection,  rad 


y  - 


a  angle  of  attack,  rad 
q  pitch  rate,  rad/ sec 


We  should  point  out  that  the  airframe  is  unstable  with  an  unstable  pole  at  s  =  +2.372. 


4.2.2  Actuator  Model 

For  preliminary  design  a  simple  actuator  model  is  desired,  but  it  should  still  be 
reasonably  accurate.  For  this  example,  the  actuator  model  represented  in  state  space  form 
is  as  follows: 


z  —  A22  +  Bjr 


(4.18) 


U>  =  C22  +  D2I: 


(4.19) 


where 


A2 


-13  0 

0  -13 


(4.20) 
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b2 


Cj  = 


D,= 


13  0 

0  13 

1  0 
0  1 

0  0 
0  0 


(4.21) 


(4.22) 


(4.23) 


and 


~2 


actuator  state 


actuator  state 


v  = 


S', 


commanded  horizontal  tail  deflection,  rad 
commanded  flap  deflection,  rad 


7r  — 


Sr 

Si 


horizontal  tail  deflection,  rad 
flap  deflection,  rad 


For  this  example  this  actuator  model  is  used  to  relate  the  actual  surface  deflection  vector, 

f-n  f  /  T 

[<5e,6f] 1 ,  to  the  commanded  surface  deflection  vector,  [6e,6f]  • 
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Airframe— Plus— Actuator  Model 


A*.  ®x] 


1  V 


,  B.] 


c,  D, 


Figure  4.22:  The  airframe-plus-actuator  model 


The  actuator  model  is  inserted  in  series  with  the  airframe  model,  as  shown  in  Figure 

4.22,  calling  it  the  airframe-plus— actuator  model,  which  is  the  plant  where  v  is  the 

two— input  vector  and  y  is  the  two-output  vector.  Denote  the  plant,  the  airframe- 

FA  Bl 

plus— actuator  model,  by  the  state  space  representation  as  S  =  q  ^  ,  where 


A2  0 


BiC2  Ai 


c  =  DiCj  C,  D  =  DiDj 


(4.24) 


Thus  the  plant  in  transfer  function  matrix  form  is  given  by 


P(s)  =  C[»I-AJ-*B+D 


(4.25) 


where 


Ph(j)  Pi»(*) 
P21 (»)  P«(*) 

M 


(4.26) 
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pu(s)  =  1.9898a3  -f  276.03a2  +  4.543a  +  0.82539 

pi2(a)  =  2.1378a3  +  66.103a2  +  1.0805a  +  0.19346 

pjl(3)  =  274.25a3  +  303.54a2  +  7.6929a 

p„(a)  =  63.797a3  +  90.402a2  +  2.0213a 

p(a)  =  a5  +  15.138a4  +  17.174a3  -  138.3a2  -  2.1981a  -  0.44993 

/  /  T 

Here  the  input  is  [6e,6 fj  and  the  output  is  [a,q]  .  We  would  like  to  point  out  that  the 
plant  checks  out  to  be  both  controllable  and  observable. 

4.2.4  Stabilisation  of  the  plant 


Figure  4.23:  Block  diagram  of  the  stabilised  plant 


Figure  4.24:  The  control  system  design  of  YF— 16  CCV 
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Since  the  proposed  technique  in  this  report  requires  the  plant  to  be  stable,  it  is  necessary  to 

stabilize  the  airframe-plus-actuator  model  with  an  inner  loop.  Here  we  choose  to  feed 

back  sensed  pitch  rate,  q  only,  because  the  measurement  of  the  pitch  rate  is  much  more 

accurate  than  of  the  angle  of  attack.  This  introduces  an  additional  control  element  in  the 

design.  The  inner  loop  control  element  was  placed  in  the  forward  path  in  series  with  the 

airframe-plus-actuator  as  shown  in  Figure  4.23.  By  a  careful  analysis  of  the  pole-zero 

pattern  of  the  plant,  a  stabilizing  controller  T(s)  was  designed,  which  is  given  by 

_  „  5a2  +  0.0813a  +  0.11283 

T(s)  = - - - 

a  +  0.02a  (4.27) 

This  controller  is  shown  in  the  block  diagram  in  Figure  4.23.  After  doing  some  routine 
block  diagram  manipulation,  the  resultant  stabilized  plant  comes  out  to  be 

0n(s)  9u(s) 

92\{s)  9h(s) 

G(5)  =  (4.28) 

where 

5U(a)  =  9.9489a6  +  1509.6a5  +  17989a4  +  625.59a3  +  463.91a2  +  7.6291a  +  1.2107 

gl3(s)  -  2.1379a8  +  93.937a5  +  3159.1a4  +  94.808a3  +  55.052a2  +  0.63732s 

521(3)  =  1371.2a6  +  19366a5  +  20114a4  +  1258a3  +  454.23a2  +  11.284s 

522(a)  =  63.797a®  921.05a5  -f  1195.6a4  4-  49.822a3  +  0.52553a2 

5(a)  =  a8  +  28.158aT  +  1585.8a6  +  19455a6  +  18316a4  +  1193a3  +  447.8a2  +  11.167a 

rr  T 

Here  the  input  is  [u i,U2j  and  the  output  is  [a,q]  .  Now,  with  the  stabilized  plant,  G(s) 
and  the  feedback  element  taken  as  H(s)  =  I2,  we  can  use  the  proposed  design  technique  to 
design  a  controller  C(s)  as  shown  in  Figure  4.24. 
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4.2.5  Design  Requirements 

Before  we  use  the  proposed  design  technique,  it  is  necessary  to  specify  the  ‘desired’ 
closed  loop  transfer  function  matrix  D(s).  The  desired  dynamics  for  the  angle  of  attack,  a, 
and  pitch  rate,  q,  to  pilot  command,  c,  for  the  flight  control  system  are  given  as  follows: 


a  21.2 


c  a3  +  7a  +  25 

(4.29) 

q  _  21.1(3  +  1.08) 
c  a3  +  7s  +  25 

(4.30) 

All  the  specifications  for  different  cases  for  this  example  were  provided  to  us  by  Mr.  Tom 
Gentry,  the  Project  Engineer.  We  wish  to  give  some  explanation  of  the  above  desired 
dynamics  based  on  the  flying  qualities  of  this  aircraft  which  are  well  documented  in  [4]  and 
our  description  below  is  taken  from  the  above  paper. 

MIL-F-8785C  [21]  is  used  to  determine  the  basic  dynamic  requirements.  We 
concentrate  here  on  the  short-period  and  time  delay  requirements  to  achieve  desired 
performance  (Level  1).  For  a  given  aircraft  Class  (here  Class  IV,  fighter)  and  Flight  Phase 
Category  (here  Category  A,  precise  control  and  gross  maneuvering  up— and— away)  the 
appropriate  equivalent  short— period  requirements  (Sections  3.2.2. 1 . 1  and  3.2.2. 1.2)  and 
equivalent  time  delay  requirement  (Section  3.5.3)  for  Level  1  flying  qualities  are 
determined.  These  equivalent-system  parameters  are  to  be  obtained  by  matching  classical 
responses  to  the  actual  pitch— rate  and  normal-acceleration  transfer  functions.  The  section 
on  equivalent  short— period  natural  frequency  also  requires  the  use  of  the  parameter  n fa  in 
g’s  per  rad.  This  parameter  is  the  ratio  of  the  change  in  a  steady-state  normal  load  factor 
(— az/g)  per  unit  change  in  angle  of  attack.  These  are  dynamic  requirements  which  we  will 
use  to  define  the  dynamic  responses  of  the  design. 

At  0.8  Mach  and  20,000  feet  altitude,  air-to-air  combat  is  an  appropriate  task. 
This  task  is  in  Flight  Phase  Category  A  in  MIL-F— 8785C.  The  factor  n/a  may  be 
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estimated  from  the  larger  of  the  two  non— zero  numerator  rate  to  elevator  deflection 
transfer  function.  Here  this  factor  is  1.082.  1/T^,  is  approximately  equal  to  L^,  the 
dimensional  derivative.  From  [8],  for  example, 

1  «  V  1  (4.31) 

where  V  is  the  aircraft  speed  in  feet  per  second  and  g  is  the  acceleration  due  to  gravity  in 
ft/sec2.  Here  V  is  829.6  ft/sec  and  g  is  32.17  ft/sec2.  The  parameter  n/a  is  estimated  to 
be  27.9  g’s/rad. 

Based  on  this  value  of  n/a,  Figure  1  in  MIL— F— 8785C  for  Flight  Phase  Category  A 

shows  an  allowable  range  for  Level  1  flying  qualities  of  0.28  and  3.6  for  the  parameter 
2 

u)s p/(n/a).  This  parameter  is  referred  to  as  the  Control  Anticipation  Parameter  (CAP). 
For  n/a  =  27.9  and  the  allowable  range  of  CAP,  an  allowable  range  of  equivalent  short- 
period  natural  frequency,  usp,  is  2.79  to  10  rad/sec.  A  design  value  of  5  rad/sec  was  chosen 
here  as  being  a  conservative  choice  to  allow  for  changes  up  or  down  and  still  be  in 
compliance  with  the  specification. 

The  allowable  range  of  equivalent  short— period  damping  ratio  is  found  in  Table  IV 
in  MIL— F-8785C.  For  Flight  Phase  Category  A,  Level  1  flying  qualities,  this  table 
permits  a  damping  ratio  of  0.35  to  1.30.  A  value  of  0.7  is  chosen  as  a  reasonable  value  for 
the  design.  This  value  is  also  typical  of  many  real  aircraft. 

The  last  parameter  to  be  specified  is  equivalent  time  delay.  For  Level  1  flying 
qualities  in  all  Flight  Phase  Categories,  the  maximum  equivalent  time  delay  allowed  is  0.1 
sec.  in  Table  XIV  [21].  Equivalent  time  delay  is  the  total  of  pure  time  delay  plus  the 
contributions  of  high-frequency  elements  such  as  actuators,  compensation  filters,  etc.  in 
the  equivalent  system  match  of  the  aircraft’s  actual  dynamics.  The  intent  of  this 
requirement  is  to  minimize  equivalent  time  delay,  so  for  our  model  response  we  will  take  a 
value  of  0  (zero). 
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Thus,  based  on  the  task  requirement  and  the  aircraft  class,  the  design  parameters 
are  an  equivalent  short— period  natural  frequency  of  5  rad/sec,  a  damping  ratio  of  0.7  and 
no  equivalent  time  delay. 

There  is  now  one  last  detail  of  the  design  requirements  to  specify,  which  is  to 
completely  define  the  exact  dynamic  specification  for  the  aircraft  response  of  interest  for 
the  design  process  to  follow  This  is  not  always  as  straightforward  as  it  may  seem.  The 
two  model  parameters  previously  specified,  short-period  natural  frequency  and  damping 
ratio,  only  define  the  denominator  dynamics  of  the  short— period  approximation  to  the 
aircraft  dynamics.  The  numerator  dynamics  must  be  specified  also.  For  short— period 
dynamics,  this  requires  that  the  pitch  rate  numerator,  for  example,  be  first-order,  which 
for  a  conventional  aircraft  is  the  factor  at  1/T^.  As  we  have  noted,  1/T^  is  largely 
determined  by  L^,  which  in  turn  is  determined  by  the  choice  of  wing  size  and  planform. 
Here  we  assume  that  this  choice  has  been  made  to  satisfy  other  design  requirements  and 
will  check  its  flying  qualities  implications. 

All  of  the  longitudinal  responses  are  coupled  through  the  airframe.  Because  only 
one  control  surface  is  available  to  control  the  longitudinal  responses  of  the  aircraft, 
complete  specification  of  one  response,  such  as  pitch  rate  or  normal  acceleration, 
completely  specifies  all  other  responses  as  coupled  by  the  particular  airframe. 

As  specified,  these  desired  dynamics  are  primarily  generated  by  the  elevator 
deflection,  6e  In  view  of  this,  we  decided  to  place  the  desired  dynamics  in  the  first  column 
of  D'(s)  to  indicate  the  major  elevator  activity.  Meanwhile  the  second  column  of  D'(s) 
should  have  insignificant  responses.  Therefore,  we  let  a/c  and  q/c  be  the  (1,1)  and  (2,1) 
element  of  D'(s),  respectively.  Since  we  want  the  responses  of  the  second  column  elements 
of  D'(s)  to  be  insignificant  and  keeping  in  view  the  degree  constraints  for  ‘good  fit’,  we 
take  the  transfer  function  for  each  of  the  (1,2)  and  (2,2)  elements  in  D'(s)  as  l/(s+50)2. 
Thus, 
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D'(*)  = 


s-  +  7a  +  25  (a  +  50)2 
21.1(3  +  1.08)  l 

J2  +  7s  +  25  (7+  50)* 


(4.32) 


In  order  to  satisfy  the  degree  condition  as  in  remark  5  in  Chapter  3,  we  enhance  D'(s)  by 
multiplying  the  (2,1)  element  by  50/(s+50).  Thus,  the  ‘desired’  closed  loop  transfer 
function  matrix  is  taken  as 


21. 2d2  +  2120a  +  53000  a2  +  7j  +  25 

1055a2  +  53889.4a  +  56970  a2  +  7a  +  25 
=  a4  +  107a3  +  3225a2  +  20000a  +  62500 


(4.33) 


4.2.6  The  Design 

In  this  example,  the  matching  of  the  first  column  elements  of  D(s)  is  of  paramount 
importance.  The  bandwidth  of  the  (1,1)  element  is  found  to  be  [0,5]  rad/sec,  and  of  the 
(2,1)  element  is  found  to  be  [2.5,8]  rad/sec.  Therefore  the  frequency  interval  of  [0,8]  or  less 
should  be  considered  for  the  minimization  of  the  error  function  Em.  Using  the  program 
with  the  plant  in  Equation  (4.28)  and  the  ‘desired’  transfer  function  matrix  in  Equation 
(4.33),  we  varied  the  frequency  interval  and  the  order  of  the  controller  until  we  obtained 
satisfactory  results.  Fairly  satisfactory  results  were  obtained  by  using  the  frequency 
interval,  [0,2]  for  1  <  i,  k  <  2,  with  a  second  order  controller.  The  proposed  technique 
yielded  the  following  transfer  function: 


C(a)  = 


cn(a)  cu(a) 
c2i(a)  c„(a) 


(4.34) 
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where 

c„(s)  =  -0.14952j2  +  19.85508s  +  19. 15416 
=  -0.14952(s  +  0.95779)(s  -  133.75) 
cia(j)  =  -0.0010493s2  +  0.0099626s  +  0.008956 
=  -0.0010493(s  +  0.82694)(s  -  10.321) 
cai(s)  =  2.23717s2  +  1.31541s  +  1.79723 
=  2.23717 (sJ  +  0.58798s  +  0.80335) 
cjj(s)  =  0.022089s2  +  0.0003915s  +  0.003185 
0.022089(s2  +  0.017727s  +  0.1442) 
c(s)  —  s2  -f  6.6078s  +  3.52986 
-  (s  +  0.5862)(s  +  6.0216) 

Note  that  the  controller  C(s)  is  stable.  Also  the  compensated  closed  loop  transfer  function 
matrix  is  given  by 

F(s)  =  [I2  +  G(s)C(s)Hrs)]-,G(s)C(sl  (4.35) 

fn(s)  Ms) 
fn{s)  fn{s) 

. 7(a)  (4  36) 

where 

/u(s)  T  3.295 17(s2  |  0.016392s  +  0.0029554)(s2  4  0.016302s  1  O.02247) 

x (s  +  0.58597)(s  +  0.98225)(s  +  6.0649)(s  +  10.728)(s2  4  44.29s  +  9988.9 
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fu(fi)  -  0.036783(s2  I  0.0163C5.<i  \  0.0028624)(.s2  +  0.016619a  -  0.022267) 

x(a  +  0.5862)(s2  +  2.9794a  +  2.6362)(a  +  6.0216)(a2  +  13.019a  +  1754.3) 
/31(a)  — 62.302a(a  +  0.025953)(a2  +  0.016325a  +  0.022445)(a  +  0.5862) 

x(a  +  0.98307)(a  +  1.059)(a  +  6.0216)(a  +  13)(a  -  440.03) 

/22(a)  =  — 0.029687a(a  +  0.029317)(a2  +  0.017013a  +  0.022136)(a  -f  0.29824) 
x(a  +  0.97521)(aJ  +  1.4627a  +  6.5752)(a  +  18.807)(a  -  476.78) 

/(a)  =  (a2  +  0.017205a  +  0.002466l)(a2  +  0.016417a  +  0.022829)(a  +  0.58597) 
x  (a  +  0.98263)(a2  +  7.0979a  +  24.672)(a  +  6.0645)(a  +  12.471) 
x(a2  +  14.138a  +  1353.2) 


4.2.7  Design  Comparison 

The  frequency  responses  of  each  element  of  the  compensated  and  the  ‘desired’  closed 
loop  transfer  function  matrices  are  given  in  Figures  4.25  through  4.32.  Since  the  magnitude 
responses  of  the  compensated  and  the  ‘desired’  system  of  the  (1,2)  and  (2,2)  elements  are 
relatively  very  small  as  shown  in  Figure  4.26  and  Figure  4.28,  therefore  the  matching  of  the 
phase  responses  of  the  elements  are  not  of  interest.  And  also,  the  time  responses  due  to  the 
unit,  step  inputs  j^Q  j  of  the  compensated  and  the  ‘desired’  systems  are  given  in  Figures 
4.33  and  4.34.  Now  the  transfer  function  matrix  of  the  elevator  and  flap  deflections  of  the 
airframe  with  respect  to  the  pilot  commanded  input  is  given  by 


where 


L(-) 


ln(a)  Ini* 

hiif)  J 


l(s) 


(4.37) 
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/u(a)  =  —  9.7l88(a2  4  0.016522a  4  0.0030627)(a2  -  0.01626a  -  0.022566) 
x(a  4  0.58561)(a  +  0.98206)(a  -  3.1956;  a  4  5.7109)(a  -  126  39) 
x(a2  +  11.846a  4  49.841) 

l17{s)  =  —0.068207(aJ  4  0.016975a  +  0.0027477j(a2  4  0.01626  +  0.022566) 
x(a  4  0.93667)(a2  +  0.21863a  4  0.26736)(a2  4  7.4103a  4  26.911) 
x(a2  4-  3.6697a  4  1179.2) 

/21(a)  =  29.0836a(a  +  0.026042)(a2  +  0.014474a  +  0.024529)(a  +  0.57672) 
x(a  4  0.97889)(a2  4  0.59512a  4  0.83605)(a  4  5.S277) 
x (a2  4  14.335a  4  1376.4) 

f„(a)  =  0.2872(a2  +  0.01738a  -  0.0023487)(a2  4  0.016484a  4  0.023031) 

x(a  4  0.9656 ")(a2  -  0.00002985a  4  0.11558)(a2  +  6.6353a  4  23.155) 
x(a2  4  14.149a  4  1354.5) 

l(s)  =  (a2  4  0.017205a  4  0.0024661)(a2  +  0.016417a  4  0.022829)(a  4  0.58597) 
x(a  4  0.98253)(a2  -+•  7.0979a  4  24.672)(a  4  6.0645)(a  4  12.471) 
x (a2  4  14.138a  4  1353.2) 

The  time  response  of  the  elevator  and  flap  deflections  of  the  airframe  due  to  the  unit  step 
command  input  J  of  the  compensated  system  are  also  given  in  Figures  4.35  and  4.36. 
These  figures,  indeed,  bring  out  the  expected  results  that  the  elevator  deflection  is  the 
major  control  surface  for  the  desired  dynamics  of  the  aircraft. 

Comparison  of  the  magnitude,  phase  and  time  responses  of  the  compensated  system 
with  C(s)  with  the  ‘desired’  system  as  shown  in  Figures  4.25  through  4.34  shows  that  the 
results  obtained  from  the  proposed  design  technique  are  excellent. 
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4.3  YF-16  CCV.  Case  H 


We  pursue  the  same  aircraft  example  viz.  YF-16  CCV  whose  model  and  actuator 
for  longitudinal  control  are  given  in  (4.2).  In  this  case,  Mr.  Tom  Gentry,  the  Project 
Engineer,  gave  us  the  following  ‘desired’  dynamics: 


a  _  13. 1  rad 

c  -  s2+7.7s+30.25  unit 


(4.38) 


q_  13.0(s+1.25)  rad/sec 
c  ~  s2+7. ys+30.25  unit 


(4.39) 


The  justification  for  these  desired  dynamics  is  similar  to  the  one  in  (4.2)  and  we 
don’t  want  to  dwell  on  that.  For  reasons  similar  to  the  ones  for  (4.2),  we  keep  the  desired 
dynamics  in  the  first  column  of  the  desired  transfer  function  matrix.  Keeping  in  view  the 
degree  condition  (remark  5,  Chapter  3)  as  elaborated  in  (4.2),  we  select  the  following 
transfer  function  matrix  D(s)  as  the  ‘desired’  closed  loop  system: 


D(s)  = 


13  . 1  1 

s 2+7. 7s +30. 25  (s+50 ) 2 

13.0(s+1.25)  50  1 

.  s2+TTslilT.25  •  s+SO  (s+50) 2 


(4.40) 


After  a  few  ‘trial  and  error’  attempts,  we  settled  on  a  second  order  controller  with  [0,4]  as 
the  interval  of  integration  for  all  components  in  the  error  function.  The  proposed 
technique  yielded  the  following  transfer  function  matrix  for  the  controller  C(s): 


where 


(4.41) 
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Cn(s)  =  0.0020384s2  +  0.69268s  +  0.97256 

=  0.0020384  (s+1.4099)(s+338.41) 

Ci2(s)  =  -0.000026447s2  +  0.00039187s  +  0.0012705 

=  -0.000026447(s+2.7367)(s— 17.554) 

c2i(s)  =  0.07805s2 -0.0032581s -0.87825  (4.42) 

=  0.078085  (s+3.1516)(s— 3.5688) 
c22(s)  =  -0.00033998s2  +  0.00012675s  -  0.025865 

-0.00033998  (s— 0.18641  ±  8.7203j) 
c(s)  =  0.061069s2  +  0.41695s  +  1 

=  0.061069  (s+3.4137  ±  2.1729j) 

Both  C(s)  and  the  closed  loop  transfer  function  matrix  F(s)  are  stable  for  this  design. 

The  frequency  responses  of  each  element  of  the  compensated  and  the  ‘desired’  closed 
loop  transfer  function  matrices  are  given  in  figures  4.37  through  4.44.  Since  the  magnitude 
responses  of  the  (1,2)  and  (2,2)  elements  in  the  ‘desired’  transfer  function  matrix  were 
almost  zero  and  the  corresponding  responses  in  the  compensated  system  turn  out  to  be 
small  (Figures  4.38  and  4.40),  no  attention  or  significance  is  attached  to  the  phase- 
mismatch  for  the  (1,2)  and  (2,2)  elements.  However,  the  match  for  the  magnitude  and 
phase  responses  for  (1,1)  and  (2,1)  elements  (Figures  4.33,  4.39,  4.41,  4.43)  is  excellent 
showing  thereby  the  second  order  controller  obtained  by  the  proposed  technique  gives 
excellent  results.  The  time  responses  due  to  unit  step  input  ^  j  of  the  compensated  and 
the  ‘desired’  systems  are  shown  in  Figures  4.45  and  4.46.  Again  good  results  are  confirmed 
by  the  time  response. 

4.4  YF-16  CCV.  Case  HI 

For  the  same  aircraft  model  studied  in  (4.2)  and  (4.3),  we  design  a  control  system 
by  using  the  proposed  itUinique  to  satisfy  another  set  of  specifications  given  to  us  by  Mr. 
Tom  Gentry,  the  Project  Engineer. 
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The  desired  dynamics  are: 


q  _  1 

c  ~  s2+7s+25 


(4.43) 


(4.44) 


where  7  is  the  flight  path  angle.  With  q  and  7  as  outputs,  and  <5e,  the  horizontal  tail 
deflection  and  ft,  the  flap  deflection  as  inputs,  the  transfer  function  matrix  of  the  plant  for 
the  aircraft  of  (4.2)  is  derived  as 


where 


p(s)  = 


Pus  P  i2(s 
P 2l(S  )  P22(S 


(4.45) 


Pu(s)  =  274.25s3  4-  303.54s2  +  7.6929s 

Pi2(s)  =  63.797s3  +  90.402s2  +  2.0213s 

P2i(s)  =  1.9898s3  +  1.7819s2  -  299s  -  6.8676  (4.46) 

P22(s)  =  2.1378s3  +  2.3051s2  -  89.322s  -  1.8278 

P(s)  =  s*  +  15.138s4  +  17.174s3  -  138.3s2  -  2.1981s  -  0.4493 


We,  first,  design  a  stabilizing  controller  T(s)  as  shown  in  Figure  4.47  and  given  by 


rr/^_5(s2  +  0.01626s  +  0.022566) 
T(S'  “  II  s  TO  702) - 


(4.47) 


“2. 


9 


u2 


)  > 

T(s) 

P(s) 

JL- 

• - J- 

Figure  4.47:  Stabilisation  of  the  Plant 
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With  T(s)  in  the  loop  as  in  Figure  4.47  the  closed  loop  transfer  function  matrix  G(x)  is 
calculated  as 

(4.48) 

where 

qu(s)  -  1371.2s6  +  19366s5  +  20114s<  +  1258s5  +  454. 23s5  +  11.284s 
gi2(s)  =  63.797s6  +  921.05s5  +  1195.6s<  +  49.822s5  +  ,052553s5 

g2i(s)  as  9.9489s6  +  138.41s5  -  1376.7s<  -  19488s5  -  794.08s5  -  446.6s  -  10.073  (4.49) 

g22(s)  =  2.1378s6  +  30.14s5  +  2238s<  -  1100.8s5  +  5.2304s5  +  0.11178s 
g(s)  =  s«  +  28.158s7  +  1585.8s6  +  19455s5  +  18316s4  +  1193s5  +  447.8s5  +  11.167s 
Keeping  the  degree  condition  (remark  5,  chapter  3)  in  mind,  the  ‘desired’  transfer  function 
matrix  is  taken  as 


By  using  the  proposed  technique,  a  second  order  controller  was  designed  such  that 
the  controller  and  the  compensated  closed  loop  system  are  stable  and  the  frequency 
responses  of  the  compensated  system  and  the  ‘desired’  system  are  close.  The  controller 
transfer  function  matrix  C(s)  by  the  proposed  technique  is  given  by 


where 

Cn(s) 

Ci2(s) 

c2i(s) 

c(s) 


=  0.000072834  (s  +  5.1914  ±  23.592j) 

=  -0.000035781  (s  +  0.14002  ±  9.2895j) 
-  -0.00096978  (s  +  3.4304  ±  5.8456j) 

=  0.040639  (s  +  3.5046  ±  3.5106j) 


(4.51) 


(4.52) 
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The  magnitude  and  phase  responses  of  the  compensated  and  the  ‘desired’  transfer 
function  matrices  are  shown  in  Figures  4.48  through  4.55.  Except  for  not  so  good  a  fit  at 
very  low  frequency  in  Figure  4.48,  these  plots  clearly  demonstrate  that  the  design  carried 
out  by  the  proposed  technique  gives  good  results. 

Summary 

Two  numerical  illustrative  examples  have  been  studied  in  this  chapter  which  show 
that  the  design  carried  out  by  the  proposed  technique  gives  excellent  results. 
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5.  ROBUSTNESS  ANALYSIS  FOR  YF-16  CCV 


In  this  chapter,  we  carry  out  a  modest  robustness  analysis  for  the  aircraft  of 
example  2.  The  best  way  to  approach  this  would  be  to  vary  several  parameters  in  the 
model  of  the  plant  and  then  analyze  the  stability  and  performance  robustness  of  the  design. 
But  this  is  quite  involved  and  quite  time  consuming.  We  thus  settled  for  a  less  lofty  goal. 
We  decided  to  vary  one  parameter  in  the  characteristic  polynomial  of  the  plant  after  it  was 
stabilized. 

The  marginally  stabilized  plant  for  the  F— 16  CCV  for  example  is  given  by  Equation 
(4.28).  The  characteristic  polynomial  g(s)  is  given  by 

g(s)  =  s8  +  28.158s7  +  1585.8s8  +  19455s8  +  183l6s«  +  1193s8  +  447.8s7  +  11.167s  (5.1) 

We  decided  to  study  the  stability  and  performance  robustness  of  the  compensated 
system  by  varying  the  coefficient  of  s7  in  (5.1).  We  denote  this  coefficient  by  6.  Thus  the 
nominal  value  of  6  is  28.158.  We  found  that  for  6  in  the  interval  [20,34],  the  inner  loop  of 
the  control  system  in  Figure  4.3  stays  marginally  stable.  With  the  controller  C(s)  designed 
for  the  nominal  plant  in  Example  2  and  given  by  Equation  (4.34),  the  overall  closed  loop 
system  represented  by  Figure  4.24  also  remains  stable.  Thus  for  6  in  [20,34]  both  the  inner 
loop  and  the  overall  system  remain  stable,  thereby  establishing  the  stability  robustness  of 
the  design  over  this  range  of  6.  These  results  have  been  established  by  calculating  the  poles 
of  the  system  by  varying  <5  over  [20,34]  with  a  fairly  small  step  size. 

Now  the  performance  robustness.  Figure  5.1  shows  that  the  graph  of  the  error 
function  between  the  frequency  responses  of  the  compensated  closed  loop  system  and  the 
desired  closed  loop  system  over  the  frequency  interval  [0,10]  rad/sec  when  6  is  varied  over 
the  interval  [20,34].  The  value  of  the  error  for  the  nominal  value  of  <5(=28. 158)  is  0.375. 
From  Figure  5.1  it  is  clear  that  the  error  function  for  other  values  of  6  in  [20,34]  remains 
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between  0  and  1.  This  establishes  performance  robustness  of  the  design  when  only  6  is 
varied  over  [20,34].  To  reinforce  the  performance  robustness  of  the  design,  we  have  drawn 
the  magnitude  and  the  phase  responses  of  the  compensated  and  the  ‘desired’  transfer 
function  matrices  for  some  sample  values  of  6  in  [20,34].  In  fact,  figures  are  drawn  for  S  = 
20,  25,  30,  and  34.  These  figures  also  portray  a  good  performance  robustness  of  the  design 
when  S  is  varied.  For  these  values  of  6,  the  magnitude  and  phase  responses  are  shown  in 
Figures  5.2  through  5.33.  Unit  step  responses  are  shown  in  Figures  5.34  through  5.41. 

As  was  mentioned  earlier,  this  stability  and  performance  robustness  analysis  is  only 
a  small  part  of  the  overall  robustness  analysis  of  the  system.  But  the  results  given  by  the 
analysis  given  in  this  report  are  fairly  encouraging  and  make  one  optimistic  about  the 
robustness  of  the  design  by  the  proposed  technique. 

A  complete  robustness  analysis  of  the  proposed  technique  is  left  for  future  research 

work. 

Summary 

In  this  chapter,  a  minor  robustness  analysis  for  stability  and  performance  for  the 
aircraft  example  has  been  carried  out.  It  is  observed  that  the  design  carried  out  by  the 
frequency  matching  technique  of  this  project  gives  fairly  decent  robustness  results  when 
only  one  parameter  in  the  plant  is  varied. 
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Magnitude 


Figure  4.1:  Linear  magnitude  response  of  the  (1,1)  element  of  the  closed  loop  transfer 
function  matrix  with  zero  order  controller  for  Example  1 
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Figure  4.3:  Linear  magnitude  response  of  the  (2,1)  element  of  the  closed  loop  transfer 
function  matrix  with  zero  order  controller  for  Example  1 


73 


Magnitude 


Figure  4.4:  Linear  magnitude  response  of  the  (2,2)  element  of  the  closed  loop  transfer 
function  matrix  with  aero  order  controller  for  Example  1 
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Figure  4.5:  Phase  response  of  the  (1,1)  element  of  the  closed  loop  transfer  function 
matrix  with  zero  order  controller  for  Example  1 
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Figure  4.7:  Phase  response  of  the  (2,1)  element  of  the  closed  loop  transfer  function 
matrix  with  aero  order  controller  for  Example  1 
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Figure  4.8:  Phase  response  of  the  (2,2)  element  of  the  closed  loop  transfer  function 
matrix  with  aero  order  controller  for  Example  1 
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Output 


Figure  4.9:  Time  response  of  the  first  output  of  the  closed  loop  system  with  sero 
order  controller  for  Example  1 
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Output 


Figure  4.10:  Time  response  of  the  second  output  of  the  closed  loop  system  with  tero 
order  controller  fo'  Example  1 
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Figure  4.11:  Linear  magnitude  response  of  the  (1,1)  element  of  the  closed  loop  transfer 
function  matrix  with  first  order  controller  for  Example  1 
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Figure  4.12:  Linear  rnagnjt’ide  response  of  the  (1,2)  element  of  the  closed  loop  transfer 
function  matrix  with  first  order  controller  for  Example  1 
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Figure  4.13:  Linear  magnitude  response  of  the  (2,1)  element  of  the  closed  loop  transfer 
function  matrix  with  first  order  controller  for  Example  1 
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Figure  4.14:  Linear  magnitude  response  of  the  (2,2)  dement  of  the  closed  loop  transfer 
function  matrix  with  first  order  controller  for  Example  1 
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Figure  4.15:  Phase  response  of  the  (1,1)  element  of  the  dosed  loop  transfer  fermion 
matrix  with  first  order  controller  for  Example  1 
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Figure  4.17:  Phase  response  of  the  (2,1)  element  of  the  dosed  loop  transfer  function 
matrix  with  first  order  controller  for  Example  1 


Figure  4.19:  Time  response  of  the  first  output  of  the  closed  loop  system  with  first 
order  controller  for  Example  1 
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Output 


Figure  4.20:  Time  response  of  the  second  output  of  the  dosed  loop  system  with  first 
order  controller  for  Example  1 
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Figure  4.25:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compen*ated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.26:  Linear  magnitude  response  of  tbe  (1,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.27:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
'desired'  closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.28:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  dosed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.30:  Phase  response  of  the  (1,2)  element  of  the  compensated  and  ‘desired’ 
closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.31:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired 
closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.32:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’ 
closed  loop  transfer  function  matrix  for  YF-16  CCV 
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Figure  4.33:  Time  response  of  the  ingle  of  attack  of  the  compensated  system  and  the 
‘desired’  system  Jo*  YF-16  CCV 
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Figure  4.34:  Time  response  of  the  pitch  rate  of  the  compensated  system  and  the 
‘desired’  system  for  YF-16  CCV 
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Figure  4.35:  Time  response  of  the  elevator  deflection  of  the  compensated  system  for 
YF-16  CCV 
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Figure  A  36:  Time  response  of  the  flap  deflection  of  the  compensated  system  for  YF-16 

ccv 
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Frequency  in  rad.  per  sec. 


Figure  4.3 1 .  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired'  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.38:  Linear  magnitude  response  of  the  (1,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.39:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.40:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.41:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.42:  Phase  response  of  the  (1,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.43:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2 
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Figure  4.44:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’  dosed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  2. 
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Figure  4.45  Time  response  of  the  (1,1)  element  of  the  compensated  system  and  the 
‘desired’  system  for  YF-16  CCV,  CASE  2. 
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Figure  4.46:  Time  response  of  the  (2,1)  element  of  the  compensated  system  and  the 
‘desired’  system  for  YF-16  CCV,  CASE  2. 
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Figure  4.48:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 
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Figure  4.51:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 
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Figure  4.52:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 


117 


Frequency  in  rad.  per  sec. 


Figure  4.53:  Phase  response  of  the  (1,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 
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Figure  4.54:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 
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Figure  4.55:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  CASE  3. 


Figure  5.1:  Error  versus  Delta. 
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Figure  5.2:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.3:  Linear  magnitude  response  of  the  (1,2)  dement  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.4:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.5:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 


125 


Phase  (deg.) 


Figure  5.6:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.7:  Phase  response  of  the  (1,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.8:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.10:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 


Magnitude 


.00055 

.0005 

.00045 

.0004 

.00035 

.0003 

.00025 


Figure  5.11:  Linear  magnitude  response  of  the  (1,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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5.12:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
l’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.13:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.14:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  dosed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.16:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  dosed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.17:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.18:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.19:  Linear  magnitude  response  of  the  (1,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 


139 


Magnitude 


Frequency  in  rad.  per  sec. 


Figure  5.20:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.21:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.22:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.23:  Phase  response  of  the  (1,2) 
loop  transfer  function  matrix  for  YF-16  C 


element  of  the  compensated  and  ‘desired’  closed 
CV,  DELTA  =  30. 
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Figure  5.24:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.25:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’  dosed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  30. 
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Figure  5.26:  Linear  magnitude  response  of  the  (1,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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&?d^I(Za[,“lI“',“de.  re8p0“<;  °5  ,he  (‘>2)  element  of  ‘he  compensated  and 
aesirea  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.28:  Linear  magnitude  response  of  the  (2,1)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.29:  Linear  magnitude  response  of  the  (2,2)  element  of  the  compensated  and 
‘desired’  closed  loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.30:  Phase  response  of  the  (1,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.31:  Phase  response  of  the  (1,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.32:  Phase  response  of  the  (2,1)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.33:  Phase  response  of  the  (2,2)  element  of  the  compensated  and  ‘desired’  closed 
loop  transfer  function  matrix  for  YF-16  CCV,  DELTA  =  34. 
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Figure  5.35:  Time  response  of  the  (1,2)  element  of  the  compensated  system  and  the 
‘desired’  system  for  YF-16  CCV,  DELTA  =  20. 
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Figure  5.36:  Time  response  of  the  (2,1)  element  of  the  compensated  system  and  the 
‘desired’  system  for  YF-16  CCV,  DELTA  =  25. 
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Figure  5.37:  Time  response  of  the  (2,2)  element 
‘desired’  system  for  YF-16  CCV,  DELTA  =  25. 
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